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Contributions  to  the  Theory  of  the  Properties  of 
Hydrogenated  Amorphous  Silicon 


BACKGROUND 

Knowledge  of  the  electronic  structure  of  crystalline 
and  amorphous  semiconductors  has  been  of  great  assistance  in 
the  development  of  electronic  devices.  Such  knowledge  is 
derived  from  first  principles  calculations,  the  results  of  which 
can  be  tested  against  experiments  and  subsequently  used  in 
interpreting  measurements  or  in  guiding  technological  applica¬ 
tions  . 

The  problem  of  improving  the  efficiency  of  solar  cells  by 
using  samples  of  hydrogenated  amorphous  silicon  (a-Si:H)  has 
received  a  lot  of  attention  recently.  Considerable  progress 
has  been  made  in  this  area  by  performing  both  experimental 
and  theoretical  studies  of  the  role  of  hydrogen  in  ^morphous 
silicon.  Very  recently  a  breakthrough  was  reported  in  the 
manufacturing  of  thin  films  of  a-Si:H  which  produced  solar  cells 
with  an  efficiency  greater  than  10%.  This  achievement  suggests 
that  the  a-Si  solar  cell  technology  now  presents  a  realistic 
alternative  to  the  crystalline  silicon  technology. 

2 

In  1975  Spear  and  LeComber  succeeded  in  substitutional- 
ly  doping  amorphous  Si  (a-Si),  produced  by  decomposition  of 
silane,  through  the  incorporation  of  phosphorous  and  boron 
impurities.  This  work  became  the  starting  point  of  the  current 
intensive  investigation  of  the  electronic  properties  of  a-Si. 
Paul  et  al.,  who  also  presented  similar  doping  results  in 
a-Si  produced  by  sputtering  in  an  Ar+H  plasma,  suggested  that 
the  origin  of  the  doping  effects  as^  well  as  other  good  elec¬ 
tronic  properties  of  this  material4’  is  the  passivation  of 
dangling  bonds  by  hydrogen.  It  has  been  shown  recently  that 
the  density  of  states  (DOS)  in  the  middle  of  the  gap  Cjan  b^ 
reduced  through  hydrogenation  to  values  as  low  as  5x10  cm 
eV.  However,  the  amount  of  hydrogen  in  these  Si-H  "alloys"  was 
found  to  be  as  much  as  100  times  larger  than  the  maximum  number 
of  dangling  bonds.  Therefore,  hydrogen,  in  addition  to  satu¬ 
rating  the  dangling  bonds,  also  resides  elsewhere  and  perhaps 
introduces  other  changes  in  the  electronic  structure  of  a-Si. 
Hydrogenation  not  only  eliminates  the  dangling  bond  states  from 
the  energy  gap  kut  also  widens  the  gap  as  demonstrated  by 
several  workers  "  using  different  experimental  techniques. 

Manuacript  approved  April  29,  1983. 
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Away  from  the  gap,  photoemission  measurements  revealed 
hydrogen  associated  states  well  within  the^ faience  band.  In  the 
conduction  band,  photoconductivity  data  suggest  the  forma¬ 
tion  of  Si-H  antibonding  states. 


There  have  been  several  theoretical  models  and  calcula¬ 
tions  in  order  to  gain  an  understanding  of  th^ge  experiments. 
We  mention  the  work  of  Allan  and  Joannopoulos  and  coworkers 
who  have  studied  Si-H  molecules  using  the  ^the-lattice  tech¬ 
nique;  the  calculations  of  Ching  et  al.  using  an  ortho- 
gonalized  linear-combination  of  atomic  orbitals  method  based  on 
a  continuous-ran^om-network ,  and  the  model  defect  method  of 
Divincenzo  et  al. 
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MODEL  OF  THE  DISORDER  (P1-P4) 

Our  model^  describes  hydrogenated  a-Si  by  constructing 
an  effective  lattice  whose  sites  may  have  probability  c  of  being 
vacant,  and  probability  1 -c  of  having  a  Si  atom.  In  addition, 
we  have  assumed  that  hydrogen  atoms  may  be  located  along  the 
lines  connecting  a  vacant  site  with  its  nearest-neighbors . 
Thus,  we  have  included  in  our  model,  at  random,  Si  sites, 
vacancy  sites,  and  sites  surrounded  by  one,  two,  three,  or  four 
hydrogen  atoms  which  saturate  the  Si  dangling  bonds.  Using  this 
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model  of  disorder,  we  have  used  the  coherent-potential  approxi¬ 
mation  (CPA)  to  perform  calculations  of  the  electronic  density 
of  states  (DOS).  The  results  (references  1-4  reprinted  herein) 
may  be  summarized  as  follows: 

a)  We  have  demonstrated  the  existence  of  dangling  bond  states 
in  the  gap. 

b)  We  have  shown  the  restoration  of  the  band  gap  upon  hydro¬ 
genation. 

c)  We  have  shown  that  for  a  20%  hydrogen  concentration  the  gap 

is  wider  by  0.4  eV  due  to  a  recession  of  the  valence  band. 

d)  We  have  identified  H  induced  peaks  in  the  DOS  of  the  valence 

band. 

e)  We  have  concluded  that  Si-H  antibonding  states  at  the  bottom 

of  the  conduction  band  have  wavefunctions  with  very  strong 

hydrogen  component. 

LOCAL  Si-H  BONDING  ( P5 ,  P6) 

We  have  carried  out  detailed  self-consistent  pseudopoten¬ 
tial  electronic  structure  calculations  on  the  hydrogen  saturated 
Si  vacancy,  (HSV)  (i.e.,  a  hydrogen  atom  attached  to  each  of  the 
four  Si  dangling  bonds)  to  learn  more  about  H-Si  bonding  in  a 
bulk  environment.  The  geometrical  model  consisted  of  a  crystal¬ 
line  array  of  HSV  "supercells".  The  local  DOS  confirmed  the  CPA 
results  that  valence  band  states  within  0.6  eV  of  the  gap  are 
strongly  depleted  by  hydrogenation  and  that  conduction  states 
just  above  the  gap  are  strongly  influenced  by  hydrogen.  We  have 
also  studied  relaxation  (symmetric  only)  of  the  HSV  by  calculat¬ 
ing  forces  (by  the  Hellman-Feynman  theorem)  and  total  energies. 
It  was  found  that  the  DOS  features  discussed  above  are  sensitive 
to  H-Si  bond  length  and  possibly  to  H-H  interactions.  These 
studies  suggest  that  the  stable  HSV  assumes  a  low  symmetry 
buckled  form  which  is  similar  to  the  ideal  (unrelaxed)  HSV  in 
local  H-Si  bonding  features. 

D.C.  CONDUCTIVITY  AND  OPTICAL  ABSORPTION  ( P7 ,  P8) 

In  addition  to  our  work  summarized  fbove,  we  have  perform¬ 
ed  calculations  of  the  transport  properties  of  a-Si:H.  We  have 
evaluated  the  d.c.  conductivity  by  the  Kubo-Greenwood  formula¬ 
tion  using  the  Green's  functions  generated  by  our  coherent 
potential  approximation  (CPA)  calculations.  Our  results  show 
that  the  mobility  of  electrons  is  greater  than  that  of  holes, 
apparently  due  to  the  strong  hydrogen  component  of  the  density 
of  states  just  above  the  gap. 

Finally,  we  have  been  able  to  calculate  the  a.c.  conductiv¬ 
ity  and  thereby  obtain  the  optical  absorption  coefficient.  This 
is  the  first  application  of  the  CPA  in  the  evaluation  of  the 
absorption  and  transport  properties  of  a-Si:H  using  a  realistic 
multi-band  model.  A  comparison  of  our  calculated  optical  ab¬ 
sorption  to  the  convoluted  DOS  shows  that  the  optical  gap  is 
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larger  than  the  DOS  gap.  This  difference  in  calculated  gaps  is 
due  to  remanent  wavevector  conservation  arising  from  short  range 
order  in  a-Si:H. 

Overall,  the  comparison  of  theoretical  calculations  with 
experimental  data  leads  to  the  following  important  conclusion: 
for  many  properties  of  a-Si:H  which  do  not  intimately  involve 
band  tail  states  or  gap  states,  knowledge  of  the  precise  atomic 
positions  are  not  necessary  for  their  calculation  and  comprehen¬ 
sion. 
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A  Slater-Koster  Hamiltonian  for  Si  is  constructed  using  four  (one  s  and  three  p )  orthogonal  orbitals  per  site.  This 
Hamiltonian  reproduces  reasonably  accurately  the  empirical  pseudopotential  band  structure  not  only  for  the  valence 
but  for  the  conduction  band  as  well.  It  also  determines  correctly  the  position  of  the  bound  level  of  an 
unreconstructed  vacancy  demonstrating  thus  that  the  effects  of  electronic  self-consistency  are  minor. 


I.  INTRODUCTION 

Koster  and  Slater1  proposed  in  1954  a  Green’s- 
function  method  for  the  calculation  of  the  changes 
in  the  electronic  structure  of  perfect  crystals 
caused  by  the  presence  of  localized  defects.  The 
first  numerical  implementation  of  this  method  was 
made  by  Callaway  and  Hughes2  on  silicon.  How¬ 
ever,  the  first-principles  numerical  evaluation  of 
the  matrix  elements  in  a  Wannier  representation 
made  this  approach  too  complicated. 

Recently,  due  to  the  work  of  Bernholc  et  al .'•* 
and  that  of  Baraff  and  Schluter,  significant  pro¬ 
gress  has  been  made  in  performing  realistic  de¬ 
fect  calculations.  These  authors  have  presented 
similar  formalisms  and  actual  calculations  which 
seem  to  indicate  a  level  of  accuracy  approaching 
that  of  band-structure  calculations  for  perfect 
crystals. 

The  purpose  of  the  present  work  is  to  show  that 
for  Si  an  accurate  orthogonal  basis  Slater-Koster7 
parametrization  can  be  obtained,  which  when  used 
in  the  Koster-Slater1  impurity  method  gives  re¬ 
sults  in  close  agreement  to  the  self-consistent 
calculations  of  Refs.  4,  5,  and  6.  The  orthogo¬ 
nality  of  the  basis  is  a  clear  numerical  advantage 
especially  in  attempts  to  treat  disordered  Si  by  an 
effective  Slater-Koster  (SK)  Hamiltonian.  On  the 
other  hand  we  recognize  that  the  orthogonality 
assumption,  which  implies  no  atomiclike  basis, 
makes  it  difficult  to  estimate  modifications  in  the 
matrix  elements  resulting  from  structural 
changes. 

Our  SK  Hamiltonian  produces  a  very  good  fit  to 
the  pseudopotential  band  structure  of  silicon,' 
and  reproduces  fairly  accurately  not  only  the 
band  gap  but  the  conduction  band  as  well,  unlike 
all  the  previous  attempts. 

In  Sec.  II  we  describe  our  Slater-Koster  inter¬ 
polation  scheme  and  in  Sec.  Ill  we  apply  it  to  the 
ideal  Si  vacancy  problem. 
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H.  SLATER-KOSTER  INTERPOLATION 

Since  Slater  and  Koster7  proposed  the  use  of  the 
tight-binding  method  as  an  interpolation  scheme 
there  have  been  many  attempts  to  apply  this  idea 
to  the  band  structure  of  silicon.3'15  Most  of  these 
calculations  give  a  good  description  of  the  valence 
band  but  the  band  gap  is  too  wide  and  the  conduc¬ 
tion  band  too  narrow.  For  example,  the  calcula¬ 
tion  of  Chadi  and  Cohen12  and  also  that  of  Chadi12 
gives  a  gap  of  3  eV  and  no  reasonable  representa¬ 
tion  of  the  conduction  band.  However,  Chadi1*  has 
succeeded  in  obtaining  the  correct  gap  and  con¬ 
duction  band  by  a  different  approach  in  which  he 
used  as  adjustable  parameters  the  exponential  de¬ 
cay  constants  of  Slater  orbitals  together  with  an 
empirical  pseudopotential  Hamiltonian.  We  con¬ 
sider  this  method  as  not  falling  into  Ihe  same 
category  as  SK  fits  and,  therefore,  we  will  not 
compare  with  it  except  to  say  that  our  results  are 
in  close  agreement. 

The  SK  fit  which  has  been  repeatedly  quoted  in 
the  literature  as  giving  a  realistic  tight-binding 
Hamiltonian  is  that  of  Pandey  and  Phillips.14  This 
calculation  overestimates  the  gap  to  1.4  eV  and 
gives  a  very  narrow  conduction  band.  As  stated 
by  the  authors,  their  calculations  are  not  reliable 
for  the  conduction  band.  We  will  give  a  detailed 
comparison  with  Pandey  and  Phillips  (PP)  after 
we  have  given  details  about  our  own  calculation. 

Our  SK  fit  uses  an  orthogonal  basis  set  of  s  and 
p  functions,  thus  our  nonsymmetrized  Hamiltonian 
is  an  8  x  8.  We  have  used,  as  adjustable  parame¬ 
ters,  20  three-center  interaction  integrals  which 
include  first,  second,  and  third  neighbors.  These 
parameters  were  determined  by  nonlinear  least- 
squares  fit  to  the  local  pseudopotential  results  of 
Pickett.'  This  step  was  done  after  reducing  the 
8x8  secular  equation  by  symmetry  following  the 
original  paper  of  Slater  and  Koster7  and  also  that 
of  Dresselhaus  and  Dresselhaus.'  This  block 
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*  Reference  14. 

tiiagonalizatlon  of  the  8  x  8  matrix  is  essential  for 
obtaining  a  reliable  fit.  We  have  fit  in  this  way 
all  eight  bands  of  Pickett's  calculation*  for  a  grid 
of  20  *  points  in  the  irreducible  Brillouin  zone. 

The  rms  fitting  error  was  less  than  0.25  eV  for 
the  valence  bands  which  is  slightly  better  than  the 
0.30-eV  value  quoted  by  PP.  For  the  conduction 
band  our  fitting  errors  were  0.38  eV  (5th  band), 
0.50  eV  (6th  band),  and  0.90  eV  (7th  and  8th  bands). 
Pandey  and  Phillips1*  do  not  give  rms  errors  for 
the  conduction  band;  our  own  estimates  using  their 
parameters  and  the  pseudopotential  results  of 
Pickett*  are  0.7  eV  (5th  band),  1.5  eV  (8th  band), 
and  3.5  eV  (7th  and  8th  bands).  It  is  clear  there¬ 
fore  that  our  own  SK  fit  is  of  supertor  accuracy. 
We  attribute  this  improvement  to  our  inclusion,  in 
contrast  to  PP,  of  the  third- neighbor  interac¬ 
tions  .  In  Table  I  we  list  our  SK  parameters 
fol  wing  the  not'  on  of  the  SK  paper7  and  com- 
pa.  ,L-n  wi'  ..ose  of  PP.  The  PP  parame¬ 
ters  •*  ”  been  converted  from  the  two-center  to 
the  three-center  notation  in  a  straightforward 


FIG.  1.  (a)  Energy  bands  for  SI  resulting  from  our 
Slater -Koster  Hamiltonian,  (b)  Energy  baods  for  Si 
from  the  pseudopotential  calculation  of  Pickett. 

manner.7  Figure  1(a)  shows  our  SK  energy  bands, 
and  Fig.  1(b),  for  comparison,  the  bands  from  the 
empirical  pseudopotential  of  Pickett.*  One  can 
see  that  we  have  even  reproduced  the  plane-wave 
character  of  the  conduction  bands  in  contrast  to 
the  narrow  conduction  states  given  by  Pandey  and 
Phillips.14  Our  band  gap  from  the  SK  bands  is  1 
eV  wide,  which  is  almost  exactly  the  value  from 
the  pseudopotential  results.* 

In  order  to  calculate  the  densities  of  states  (DOS) 
we  have  generated  from  the  SK  Hamiltonian  eigen¬ 
values  and  eigenvectors  for  89  k  points.  These 
results  were  then  used  in  the  tetrahedron  method.17 
The  resulting  DOS  are  shown  in  Fig.  2  where  one 
can  note  that  the  well-accepted  three-peak  struc¬ 
ture  of  the  valence  bands  is  reproduced  as  well 
as  a  fairly  accurate  value  for  the  gap  and  very 
reasonable  shape  for  the  conduction  bands.  Also, 
the  angular  momentum  decomposition  gives  the 
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FIG.  2. 
toaian. 


Total  and  angular  momentum-decomposed  densities  of  states  for  Si  derived  from  our  Slater-Koster  HamiU 


expected  strong  s  character  at  the  bottom  of  the 
valence  band  and  strong  p  character  at  the  top  of 
the  valence  band  and  the  conduction  band.  In  Fig. 

3  we  present,  calculated  in  the  same  way,  the 
DOS  from  the  Pandey  and  Phillips14  parameters. 
The  similarity  in  tne  valence  band  and  the  clear 
disagreement  in  the  conduction  band  is  evident. 

III.  IDEAL  VACANCY  IN  SILICON 

Bernholc  and  Pantelides3  have  given  a  clear  sum¬ 
mary  of  the  Koster-Slater*  theory,  and  performed 
calculations  of  the  ideal  vacancy  in  Si  based  on  the 
SK  parameters  of  PP.14  We  have  done  calculations 
along  the  same  lines  using  the  PP  parameters, 
for  the  purpose  of  checking  our  computer  codes, 
and  also  using  our  own  SK  parameters.  We  have 
reproduced  the  results  of  Bernholc  and  Pantelides3 
using  the  PP  parameters.  In  particular,  using  the 
PP  parameters  we  find  a  bound  state  of  T,  ( p - 
like)  symmetry  at  0.27  eV  above  the  top  of  the 
valence  band. 


However,  using  our  own  SK  parameters  we  find 
the  bound  state  of  T,  symmetry  to  be  at  0.75  eV 
above  the  top  of  the  valence  band,  i.e.,  much 
closer  to  the  conduction  band  rather  than  the 
valence  band.  The  bound  levels  are  determined  in 
general  as  solutions  of  the  equation 

detIK-  -£  <^a..(£)v...J=o, 

where  Caa..  (S)s(a J  (£  -tfTM  a")  is  the  unper¬ 
turbed  (crystalline)  Green’s  function’8  and  Pa..a. 
are  the  matrix  elements  of  the  perturbing  poten¬ 
tial.  In  the  present  case  Gaa,.(£)  is  diagonal 
when  a,  a”  refer  to  the  four  orbitals  of  the  same 
site  and  Va..„.  is  diagonal  and  infinite  when 
a”,  a'  refer  to  the  four  orbitals  of  the  vacant  site 
and  zero  otherwise.  Thus  the  equations  deter¬ 
mining  the  bound  levels  become 

G,(£)  =  0,  (1) 

G,(  £)  =  0  ,  (2) 
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FIG.  3.  Total  and  angular  momentum-decomposed  densities  of  states  for  SI  derived  using  the  Slater-Koster  param¬ 
eters  of  Pandey  and  Phillips. 


where  Gf,  Gs  are  the  p  and  s  diagonal  matrix  uncertainties  of  the  calculation  and  thus  can  be 

elements  of  G°,  respectively.  In  Fig.  4  we  plot  neglected.  Hence,  it  appears  to  js  that  the 

the  real  part  of  Gp  (ImGf  is  zero  in  the  gap)  vs  question  of  carrying  the  calculation  to  self-con- 

E,  showing  thus  the  solution  of  Eq.  (1)  at  E  sistency  is  of  secondary  importance.  What  is 

=  0.75  eV.  On  the  other  hand,  the  plot  of  the  real  needed  is  a  tight-binding  Hamiltonian  which  gives 

part  of  G,  versus  E  in  Fig.  5  shows  that  Eq.  (2)  the  correct  gap  and  a  good  representation  of  the 

has  no  solution,  i.e.,  no  s-like  bound  level  exists.  valence  as  well  as  the  conduction  band.  With  such 

The  important  point  in  this  section  is  that  the  a  Hamiltonian  as  a  starting  point  the  electronic 

position  of  the  bound  state  that  we  have  calculated  self-consistency  effects  seem  to  be  of  minor  im- 

(0.75  eV)  is  in  close  agreement  with  the  0.8-eV  portance.  On  the  other  hand,  the  effects  of  lattice 

value  reported  by  Bernholc  et  al.*  and  obtained  relaxation  due  to  the  vacancy  (i.e.,  the  effects  of 

from  electronically  self-consistent  calculations.  ionic  self-consistency)  are  quite  important  as 

Our  value  is  also  in  good  agreement  with  the  0.7-  suggested  recently  by  White  and  Ngai15  and  demon- 

eV  value  given  by  Baraff  and  Schluter5'’  and  also  strated  by  elaborate  first-principles  calculations 

by  an  electronically  self-consistent  approach.  by  Baraff  et  al.:a  and  Lipari  et  at.21 

This  confirms  the  suspicion  of  Bernholc  et  al.*  In  conclusion,  we  have  succeeded  in  construct  - 

that  the  position  of  the  bound  level  in  the  gap  is  ing  an  orthogonal  s,p 3  basis  Slater-Koster  Hamil- 

determined  mainly  by  the  crystalline  SK  Hamil-  tonian  for  Si  which  reproduces  reasonably  ac- 

tonian.  Indeed,  the  present  calculation  shows  that  curately  the  crystalline  band  structure  not  only  of 
the  effect  of  electronic  self-consistency  on  the  the  valence  band  but  of  the  conduction  band  as 

position  of  the  bound  level  is  within  the  numerical  well.  This  Hamiltonian  was  shown  also  to  be 
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FIG.  4.  The  real  part  of  the  p-like  Green  function 
Gt  plotted  as  a  function  of  energy.  Note  that  the  gap  is 
In  the  range  0«  £<  1.0  eV. 

capable  of  calculating  correctly  the  effects  of  un¬ 
reconstructed  vacancies.  In  forthcoming  publica¬ 
tions  we  employ  this  Hamiltonian  together  with  a 
coherent  potential  approximation  approach  to 
study  the  electronic  structure  of  amorphous  and 
hydrogenated  amorphous  Si.  In  such  a  compli¬ 
cated  calculation,  the  orthogonality  and  the  small 
number  (four  orbitals  per  site)  of  the  basis  of  the 
present  SK  scheme  are  very  important  advantages 


FIG.  5.  The  real  part  of  the  s-like  Green  functions 
C,  plotted  as  a  function  of  energy.  Note  that  the  gap  Is 
in  the  range  0=s  l.u  eV. 
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ABSTRACT 

We  have  used  the  coherent  potential  approximation  to  calculate 
the  electronic  densities  of  states  for  a  model  of  hydrogenated 
amorphous  silicon.  The  results  are  in  good  agreement  with  photo- 
emission,  optical,  photoconductivity  and  photoluminescence  data. 


Recent  experimental  work  in  a-S1:H  has  shown  that  hydrogen 
reduces  the  density  of  states  (DOS)  in  the  band  gap  by  several 

1  2  2 
orders  of  magnitude,  and  also  widens  the  optical  gap.  *  In  addi¬ 
tion,  other  experiments  have  shown  that  hydrogen  modifies  both  the 

4  5 

valence  and  conduction  bands. 

For  the  purpose  of  understanding  these  experiments  several 

theoretical  models  and  calculations  have  been  proposed.®  We  present 

here  a  brief  account  of  our  calculations^  which  we  believe  provide 
satisfactory  explanations  of  the  above  mentioned  measurements. 

Our  model  assumes  an  effective  lattice  whose  sites  may  have 
probability  c  of  being  vacant,  and  probability  1-c  of  having  a 
Si  atom.  In  addition,  we  have  assumed  that  hydrogen  atoms  may 
be  located  along  the  lines  connecting  a  vacant  site  with  Its  nearest 
neighbors.  Thus  we  have  included  In  our  model,  at  random.  Si  sites, 
vacancy  sites,  and  sites  that  have  one,  two,  three,  or  four  hydrogen 
atoms.  Using  this  model  of  disorder,  we  have  used  a  tight-binding 
form  of  the  coherent-potential  approximation  (CPA)  to  perform  de¬ 
tailed  calculations  of  the  electronic  DOS.  Since  we  have  not  allowed 
for  reconstruction,  we  believe  that  the  most  stable  configuration 
to  be  compared  with  experiment  is  that  in  which  a  vacancy  is  re¬ 
placed  by  a  complex  of  four  hydrogen  atoms.  This  configuration 
Is  shown  in  two  dimensions  In  Fig.  I. 

The  starting  point  of  this  calculation  Is  a  Slater-Koster 
(SK)  Hamiltonian  fit  to  the  pseudopotential  band  structure  of  crys- 

O 

talllne  SI.  The  basis  Is  orthonormal  with  four  orbitals  (one 
s-like  and  three  p-llke)  per  atom.  We  have  used  20  SI-SI  matrix 
elements  (three-center  Interactions)  that  Include  first,  second,  and 

ISSN:0094-243X/81/730130-06$1.50  Copyright  1981  American  Institute  of  Physics 
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Figure  1.  Two-dimensional  view  of  the  atom  configuration  used  In 
the  present  calculations. 

third  nearest  neighbors.  We  have  determined  the  H-H,  and  H-Si  matrix 

elements  from  the  S1H,  molecule  on  an  sp^  basis,  and  then  converted 
them  to  our  SK  basis. 

The  CPA  effective  Hamiltonian  H  has  off-diagonal  matrix  ele¬ 
ments  which  for  the  first  nearest  neighbors  are  a  virtual  crystal 
average  (VCA)  of  the  SI-SI  and  H-SI  values.  The  off-diagonal  second 
and  third  nearest  neighbor  Interactions  were  assumed  to  have  the 
SI-SI  values.  The  diagonal  matrix  elements  E$  and  form  a  diagonal 

4x4  matrix  (since  each  site  is  associated  with  four  orbitals): 


and  they  are  determined  from  the  CPA  condition: 

J  PJTJ  -  0  <*> 
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where  Pj  Is  the  probability  for  each  configuration  j,  and 

'j  ■  UjC-s.^r1  o) 

where  G„  is  a  4x4  site  diagonal  Green's  function  corresponding  to 
the  effective  Hamiltonian  fle,  and 

Uj  -  ^ j  (000)  -I  (4) 

where  Tj(OOO)  is  a  4x4,  in  general  non-diagonal,  matrix  which  con¬ 
tains  the  on-site  SK  parameters.  Having  determined  Zs>  Ip,  and  thus 
Ge  from  the  above  equations  we  can  calculate  the  DOS  from  the  stand¬ 
ard  expression: 


p(E)  -  -  |  ImTrGe  (5) 

Our  calculated  DOS  are  shown  in  Fig.  2  for  the  case  of  5%  vacancies 
which  are  all  saturated  by  hydrogen  resulting  In  a  hydrogen  concen¬ 
tration  of  20X.  From  Fig.  2  we  note  first  a  bandgap  Eg*1.4eV  that 

is  wider  by  0.4eV  than  the  corresponding  Eg«1.0eV  which  our  tight- 

binding  Hamiltonian  gives  for  the  non-hydrogenated  case.8  This 
widening  of  the  gap  Is  due  to  a  narrowing  of  the  third  peak  of  the 
DOS  at  the  top  of  the  valence  band  which  we  have  identified  to  be 

a  result  of  a  decrease  of  the  ppt  interaction  due  to  hydrogenation.7 
This  recession  of  the  valence  band  by  0.4eV  Is  In  excellent  agree* 
merit  with  the  photoemission  measurements  of  von  Roedem  et  al. 
Figure  2  also  shows  a  site  deronposltion  of  the  DOS.  Concentrating 
on  the  H  site  DOS,  we  have  Identified  hydrogen  Induced  peaks  at 
approximately  5.2  eV,  7.6  eV,  and  13.5  eV  below  the  Fermi  level. 

This  Is  also  In  good  agreement  with  the  photoemission  data.4 

We  have  further  demonstrated  the  widening  of  the  band  gap  by 
performing  a  joint  DOS  calculation  and  presenting  our  results  In  the 
same  manner  as  that  followed  in  the  analysis  of  the  measurements  of 

Q 

the  optical  absorption  coefficient  a.  The  calculated  square  root 

of  the  joint  DOS  is  assumed  to  be  proportional  to  the  quantity 

1/2  3 

(a£)  ,  and  has  been  normalized  to  the  experimental  value  at  E*4eV. 

A  comparison  with  the  measurements  of  Cody  et  a i^  Is  shown  in  Fig.  3. 
Although  the  experimental  graph  corresponds  to  a  smaller  hydrogen 
content  (16%)  than  the  calculated  one,  the  agreement  is  good  with  the 
only  discrepancy  being  that  the  theory  predicts  a  smaller  gap. 
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Figure  2.  Total  and  site  decomposed  densities  of  states  for  S1H 
with  x*0.20.  Note  that  the  SI  and  H  DOS  have  been  multi¬ 
plied  by  1-c  and  c  respectively  (c»0.05).  The  Fermi  level 
Is  located  In  the  middle  of  the  gap. 

Figure  4  shows  the  variation  of  the  band  gap  as  a  function 
of  the  hydrogen  content.  In  all  cases  shown  In  this  graph  we  are 
dealing  with  vacancies  fully  saturated  by  hydrogen.  We  not*  a 
monotonic  Increase  of  the  band  gap  with  Increasing  hydrogen  conceit 
t  rat  Ion.  This  Is  In  qualitative  agreement  with  the  experiments  of 
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ENERGY (eV) 

Figure  3.  Plot  of  the  square  root  of  the  joint  density  states 
versus  energy.  The  open,  circles  are  the  experimental 
values  of  the  quantity  (aE)1'*  (Ref.  3). 

2 

Moustakas  et  al  who  have  studied  the  dependence  of  the  optical 
gap  on  hydrogen  content  for  sputtered  hydrogenated  a-S1. 

We  should  mention  here  that  In  addition  to  the  hydrogen- 
saturated  vacancy  case  we  have  also  considered  configurations  where 
one,  two,  or  three  hydrogen  atoms  are  present,  leaving  three,  two,  or 
one  dangling  bonds  respectively.  The  results  demonstrate  the  appear¬ 
ance  of  dangling  bond  states  In  the  gap.^ 

Turning  now  to  the  conduction  band,  we  have  found  from  the 
site-decomposed  DOS  of  Fig.  2  that  the  bottom  of  the  conduction 
band  has  strong  hydrogen  character.  This  supports  the  notion  that 

photoluminescence  data**  suggest  the  formation  of  S1-H  antibonding 
states  at  the  bottom  of  the  conduction  band. 

We  wish  to  acknowledge  discussions  with  W.E.  Pickett,  T.D. 
Moustakas,  L.L.  Boyer,  B.M.  Klein,  and  6.0.  Cody.  We  are  also 
grateful  to  T.D.  Moustakas  for  providing  his  unpublished  data  shown 
In  Fig.  4.  This  work  was  supported  In  part  by  the  Solar  Energy 
Research  Institute  vis  an  interagency  agreement  with  the  US  Depart¬ 
ment  of  Energy. 
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Figure  4. 

Plot  of  band  gap 
versus  hydrogen 
content.  The  open 
circles  are  the  mea¬ 
surements  of  Ref.  2. 
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Using  a  coherent-potential  approximation,  we  calculate  the  effects  of  hydrogenation  on  the 
density  of  states  of  silicon  containing  a  finite  concentration  of  ideal  vacancies.  The  results 
reproduce  the  main  features  seen  in  the  measurements  including  the  restoration  and  the  widen¬ 
ing  of  the  gap  with  increasing  hydrogen  content. 


Hydrogenated  a-Si  has  recently  received  a  lot  of  at¬ 
tention  partly  because  of  the  possibilities  it  offers  as  a 
device  material.'*4  In  particular,  emphasis  has  been 
given  to  the  determination  of  the  density  of  electron¬ 
ic  states  (DOS)  by  various  experimental  tech- 
niques3-1  and  to  a  lesser  degree  by  theoretical  calcula¬ 
tions.'-10  As  a  result  of  these  efforts  it  is  generally 
believed  that  the  gap  states  in  a-Si  are  associated 
with  dangling  bonds;  hydrogen  is  attached  to  these 
dangling  bonds  thus  pushing  these  states  out  of  the 
gap.  The  model  we  consider  in  this  work  demon¬ 
strates,  for  the  first  time  in  a  quantitative  way.  the 
appearance  of  dangling  bond  states  and  (heir  passiva¬ 
tion  by  hydrogen,  and  also  shows  the  role  of  hydro¬ 
gen  in  widening  the  gap.  An  alternate  computational 
technique*  using  small  Si-H  molecules  (terminated  by 
Si  Bethe  lattices)  has  produced  qualitatively  similar 
conclusions  to  those  presented  here. 

The  starting  point 'in  our  calculation  is  a  Slater- 
Koster  Hamiltonian  H, ,  including  up  to  the  third- 
nearest-neighbor  interactions;  the  basis  is  orthonor¬ 
mal  with  four  orbitals  (one  s-like  and  three  p-like) 
per  atom."  This  Hamiltonian  accurately  reproduces" 
both  the  valence  and  the  conduction  bands  of  crystal¬ 
line  Si;  we  found  that  this  level  of  accuracy  is  essen¬ 
tial  for  our  calculation.  Randomness  is  introduced  by 
assuming  that  there  is  a  finite  concentration  c  of  ideal 
vacancies  which  means  that  the  diagonal  matrix  ele¬ 
ments  of  the  Hamiltonian,  at  a  given  site 

n  have  a  probability  c  to  be  infinite  and  a  probability 
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I  —  o  to  have  the  perfect-crystal  values.  This  type  of 
disorder  produces  the  essential  feature  of  dangling 
bonds  at  the  atoms  surrounding  a  vacancy.  Topologi¬ 
cal  disorder  is  simulated  by  changing  the  tight- 
binding  parameters  E^l  777  )  and  £,0,(777  >  from 
their  crystalline  Si  values"  to  the  values  1.39  and 
0.31  eV,  respectively,  as  a  result  of  varying  the 
dihedral  angle  in  a  range  which  spans  the  eclipsed 
and  staggered  configurations.  We  have  omitted  the 
reconstruction  and  relaxation  of  the  Si  lattice  which 
tends  to  push  the  dangling  bond  states  back  into  the 
bands  as  shown  by  explicit  calculations  in  the  one  va¬ 
cancy  case. 12,13  Thus  our  model  would  overestimate 
the  number  of  dangling  bond  stales  and  consequently 
the  hydrogen  required  for  their  passivation.  The  size 
of  this  overestimation  depends  on  how  much  recon¬ 
struction  has  taken  place,  which  in  turn  is  affected  by 
the  method  of  preparation  of  the  Si-H  Him.3  For  ex¬ 
ample  reconstruction  may  be  prohibited  when  the 
film  is  grown  in  the  presence  of  a  substantial  amount 
of  hydrogen.2  Finally  we  point  out  that  no  correla¬ 
tion  among  the  positions  of  the  vacancies  has  been 
included  in  the  present  calculation.  It  has  been 
shown  elsewhere"  that  the  main  effect  of  the  tenden¬ 
cy  to  .luster  the  vacancies  together  to  form  divacan¬ 
cies.  trivacancies,  and  internal  voids,  etc.,  is  to  re¬ 
place  <  by  a  smaller  effective  c*. 

Hydrogens  are  assumed  to  be  located  along  the 
lines  connecting  a  vacant  site  with  its  nearest  neigh¬ 
bors  Thus  up  to  four  hydrogens  can  be  accomodated 

2042 


16 


22 


RAPID  COMMUNICATIONS 


2043 


around  a  Si  vacancy.  The  H-Si  distance  is  taken  as  in 
the  SiH4  molecule.  Charge-density  contours  of  a  H 
saturated  Si-vacancy  calculation1’  suggest  the  ex¬ 
istence  of  higher  than  s  components  of  the  hydrogen 
wave  function;  and  show  that  this  multiple  /  wave 
function  can  be  approximated  by  a  single  s-wave 
function.  \SH) .  which  is  displaced  towards  Si  by 
about  20%.  Using  this  wave  function  IS,/).  the  ma¬ 
trix  elements  of  SiH4,14  the  results  of  Chadi17  that 
Si  wave  functions  decay  as  exp(  — 1.04/-),  and 
Mattheiss's  values  for  hydrogenic  matrix  elements, 

we  have  estimated  matrix  elements  y,(‘  “  1 . 6) 

analogous  to  those  defined  by  Hirabayashi19  with  one 
or  both. of  the  Si  sp 1  hybrids  replaced  by  |S«>.  When 
a  hydrogen  is  missing  the  corresponding  y\  is  taken 
as  infinite.  Knowing  the  matrix  elements  y,  for  each 
configuration,  we  can  transform19  to  the  original  basis 
|//i)  and  obtain  {ni\H\jn).  For  the  H-saturated  Si 
vacancy  the  matrix  elements  in  eV  are  £fl(000) 

— 8.72,  Ea ( 000 ) — 1.60,  — 3.05. 

TTT  >  “  1 -96.  f-f  TTT }  “ 0  64- ^ ( TTT  > 

—  0.98.  The  off-diagonal  disorder  in  our  model  Ham¬ 
iltonian  was  treated  in  the  virtual-crystal  approxima¬ 
tion  by  averaging  the  above  Si-H  matrix  elements 
with  those  of  Si."  We  then  employ  the  coherent- 
potential  approximation  (CPA)  to  obtain  the  DOS. 

The  CPA  calculates  the  DOS  from  a  periodic  effec¬ 
tive  Hamiltonian  H,  which  is  obtained  by  replacing 
(ns\H\sn)  b yl,(£)  and  (np\H\pn)  by  2,(£) 
where  p  stands  for  any  of  the  three  p  orbitals.  The 
unknown  quantities  (self-energies)  2,, 2,  are  deter¬ 
mined  from  the  self-consistency  condition 

{U{\-G,U)-')-0  .  (1) 

where  the  parentheses  denote  average  over  the  fol¬ 
lowing  17  configurations  at  a  given  lattice  site:  the 
site  occupied  by  Si  (one  configuration  with  probability 
I  —  c ).  the  site  being  vacant  with  no  H  (one  config¬ 
uration  with  probability  x0),  the  site  being  vacant 
with  one  dangling  bond  occupied  by  H  (four  config¬ 
urations  each  of  probability  X|/4),  the  site  being  va¬ 
cant  with  two  dangling  bonds  occupied  by  H  (six  con¬ 
figurations  each  of  probability  x2/6).  the  site  being 
vacant  with  three  dangling  bonds  occupied  by  H  (four 
configurations  each  of  probability  jrj/4),  and  the  site 
being  vacant  with  all  four  dangling  bonds  occupied  by 
H  (one  configuration  of  probability  jc4):  G,  is  the  di¬ 
agonal  4*4  matrix  (ni\kE  —  H,)~'\ui)  and  U  is  the 
configuration-dependent  4*4  matrix  («/ 1  (H  ~H,  )\m ). 
As  expected  from  symmetry  considerations,  the  ma¬ 
trix  (U{\~G,U)~')  was  found  to  be  diagonal  with 
the  three  p  matrix  elements  equal  to  each  other. 

Thus  Eq.  (1)  is  actually  reduced  to  two  equations 
which  determine  the  unknown  quantities  2,, 2,. 

Explicit  results  were  obtained  by  assuming  that  hy¬ 
drogen  occupies  dangling  bonds  in  an  uncorrelated 
way.  In  this  case  the  probabilities  xt,  (/  -0 . 4) 


are  given,  by 

4!  cx'iAc-x)*-1  (2) 

X,“  (4 -/)!/!  (4c)4 

where  x  is  the  ratio  of  hydrogen  atoms  to  lattice  sites. 
We  found  that  the  numerical  work  is  significantly  re¬ 
duced  if  the  following  additional  approximation  is 
employed:  Eq.  (I)  is  solved  for  the  limiting  cases 
x  —0.  x  —4c  obtaining  thus  2,1  £;0),  2,(  £;0). 
2j(£.4r),  and  2,(£;4c).  For  intermediate  values  of  x 
an  average  /-matrix  approximation  is  used  with 


2,<£;.x) 


2,(£;0)  +  — 2,(£;4c); 

4c 


(3) 


i-s.p  . 


Using  these  2's  we  obtain  the  4  x  4  matrix  G,  and 
the  configuration-dependent  4*4  matrices  U.  The  Si 
and  the  H  partial  DOS  are  given  explicitly  by  the 
equations 


Ps,*<  £) 


(I  -c) 


Tr[(  I  —  G,U^,)~'G,]  . 


(4) 


P„(£)--c—  TrUi-G.Ur'G,),  .  (5) 

IT 


where  the  symbol  (  ),  indicates  average  over  only 
those  configurations  associated  with  no  Si  at  the  site 
it.  The  total  DOS  is  of  course: 


p(£)  -ps,(£)  +pn(£)  .  (6) 

In  Fig.  I  we  plot  the  DOS  in  the  gap  where  one 
sees  that  by  increasing  the  hydrogen  content  x  the 
states  in  the  gap  are  reduced  drastically  in  agreement 
with  the  experimental  findings.  It  is  our  impression 


I'lCi.  I.  Density  of  states  in  the  gap  for  different  hydro¬ 
gen  concentrations  r.  The  Termi  level  was  found  to  he  at 
- 1)  4  cV  and  was  almost  invariant  with  x 
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PIG.  2.  Total  anti  site  decomposed  densities  of  states. 

Top  panel  refers  to  5%  vacancies  and  no  hydrogen,  and  the 
bottom  panel  shows  the  ease  with  all  dangling  bonds  ter¬ 
minated  by  hydrogen. 

that  the  present  calculation  is  the  first  one  to  show  in 
a  convincing  quantitative  manner,  this  restoration  of 
the  band  gap  upon  hydrogenation. 

In  Fig.  2  we  showihe  DOS  for  the  whole  spectrum 
including  both  the  valence  and  conduction  bands 
Figure  2(a)  shows  the  total  DOS  for  the  case  of  5% 
vacancies,  i.e  .  c  —  0.05  and  x  —0.0.  This  is  com¬ 
pared  to  Fig.  2(b),  where  all  (x  -0.20)  of  the  dan¬ 
gling  bonds  are  terminated  by  hydrogen.  The  H  site 
DOS  [in  Fig  2 ( b)  1  shows  two  peaks  at  approximately 
5.2  and  7.6  cV  below  the  Fermi  level  in  good  agree¬ 
ment  with  the  photoemission  data  of  von  Roedern 
ei  at  -  Furthermore  our  results  show  a  recession  of 
the  valence  band  by  about  0.4  eV  in  agreement  with 
the  above  experiment.5  By  examining  the  site 
decomposition  of  the  DOS  and  the  limiting  case  of  an 
isolated  hydrogen  impurity  we  have  determined  that 


TAIiLb  I.  Band  gap  £,  (in  cVI  as  a  function  of  hydrogen 
content  fnr  cases  where  all  dangling  bonds  arc  saturated 
li  e  .  v  -4c).  The  estimated  uncertainty  is  ±5%. 
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the  eigenstates  at  the  bottom  of  the  conduction  band 
of  Si-II  are  resonance  states  exhibiting  an  enhanced 
probability  to  find  the  electron  around  the  hydrogen 
atom.  These  slates  can  be  thought  of  as  antibonding 
H-Si  states  mixed  with  regular  conduction-band  Si 
states.  We  would  like  to  point  out  that  the  existence 
of  this  type  of  states  has  been  suggested  in  order  to 
interpret  transport  and  photoluminescence4  results  in 
highly  hydrogenated  samples. 

Finally,  we  wish  to  emphasize  that  the  present 
work  shows  a  substantial  increase  of  the  gap  upon 
hydrogenation  in  agreement  with  experiment.1-2  This 
is  demonstrated  in  Table  I  which  gives  the  size  of  the 
gap  as  a  function  of  hydrogen  content  Tor  cases 
where  no  dangling  bonds  are  present.  According  to 
the  present  theory  the  widening  of  the  gap  is  due  to 
the  following  two  effects:  (a)  The  Si— II  bond  being 
stronger  than  the  Si-Si  bond  tends  to  push  the  bot¬ 
tom  of  the  conduction  band  towards  higher  energies, 
thus  cancelling  the  opposite  effect  produced  by  the 
diagonal  hydrogen  matrix  elements,  (b)  The  differ¬ 
ence  which  is  a  measure  of 

the  ppir  interaction,  effectively  decreases  upon  hy¬ 
drogenation;  this  causes  a  narrowing2"  of  the  third 
peak  of  the  valence  band  (VB)  which  results  in  the 
recession  of  the  top  VII  stales. 

In  conclusion  our  calculations  for  Si-lf  reproduce 
with  no  adjustable  parameters,  the  main  features  of 
the  measurements  of  the  DOS,  including  the  widen¬ 
ing  of  the  gap. 
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We  have  used  the  coherent  potential  approximation  to  calculate  the  electronic  densities 
ot  states  for  a  model  of  hydrogenated  amorphous  silicon.  The  results  demonstrate  the 
restoration  and  widening  of  the  band  gap  with  increasing  hydrogen  content.  In  the 
valence  band,  excellent  agreement  with  photoemission  experiments  is  obtained.  In  the 
conduction  band  Si-H  antibonding  states  are  predicted  that  can  be  inferred  from  photo¬ 
conductivity  measurements. 


I.  INTRODUCTION 
Basic  experimental  facts 

The  current  intensive  investigation  of  hydro¬ 
genated  amorphous  silicon  originated  in  1975  when 
Spear  and  LeComber1  succeeded  to  substitutionally 
dope  amorphous  Si(a-Si),  produced  by  decomposi¬ 
tion  of  silane,  through  the  incorporation  of  phos¬ 
phorous  and  boron  impurities.  Soon  after,  it  be¬ 
came  apparent,  as  a  result  of  work  by  Paul  et  al. 2 
who  also  presented  similar  doping  results  in  a-Si 
produced  by  sputtering  in  an  Ar  +  H  plasma,  that 
the  origin  of  the  doping  effects  as  well  as  other 
good  electronic  properties  of  this  material3,4  is  the 
H  passivation  of  dangling  bonds.  It  has  been 
shown  recently,5  for  example,  that  the  density  of 
states  in  the  middle  of  the  gap  can  be  reduced 
through  hydrogenation  to  values  as  low  as  5x  1014 
cm-3eV~'.  However,  the  amounts  of  hydrogen  in 
these  Si-H  "alloys”  were  found  to  be  as  much  as 
100  times  larger  than  the  maximum  number  of 
dangling  bonds.  Therefore  hydrogen,  in  addition 
to  saturating  the  dangling  bonds,  introduces  other 
changes  in  the  electronic  structure  of  a-Si. 

Hydrogenation  not  only  eliminates  the  dangling- 
bond  states  from  the  energy  gap,  but  also  widens 
the  gap  as  demonstrated  by  Freeman  and  Paul,6  by 
Cody  et  al.,1  and  by  Goodman  et  el. 8  using  dif¬ 
ferent  experimental  techniques. 

Away  from  the  gap,  photoe.nission  measure¬ 
ments  by  von  Roedem  et  al .’  revealed  hydrogen 
associated  states  well  within  the  valence  band.  In 

24 


the  conduction  band  the  photoconductivity  data  of 
Moustakas  et  al. 10  suggest  the  formation  of  Si-H 
antibonding  states.  The  role  of  hydrogen  in  modi¬ 
fying  the  network  is  investigated  through  a  variety 
of  experimental  techniques,  such  as  infrared  and 
Raman  spectroscopy,1  112  nuclear  magnetic  reso¬ 
nance,13  small-angle  x-ray  scattering,14  H  implanta¬ 
tion  in  c-Si,15  and  neutron  scattering  measure¬ 
ments.16  For  review  of  experimental  work  the 
reader  is  referred  to  the  articles  by  Spear,17  Mous¬ 
takas,3  and  Fntzsche.4 


Present  physical  understanding 

Unhydrogenated  a-Si  is  thought  of  3S  a  random 
network  where  the  local  tetrahedral  arrangement, 
with  bond  lengths  almost  identical  to  those  of  the 
crystalline  state,  is  retained  to  a  high  degree.  An 
idealization  of  this  concept  is  the  so-called  “ideal 
random  network,”  where  the  tetrahedral  coordina¬ 
tion  is  satisfied  throughout  with  very  small  fluctua¬ 
tions  in  the  bond  lengths  and  larger  fluctuations  in 
other  longer  range  geometrical  aspects  such  as  ring 
sizes,  dihedral  angles,  etc.  This  ideal  random  net¬ 
work  defines  the  concept  of  topological  disorder. 

In  reality,  there  are  important  deviations  from  the 
ideal  network,  such  as  vacancies  and  other  strong 
local  distortions,  which  may  even  cluster  together 
to  form  voids,  internal  surfaces,  etc.  Associated 
with  this  type  of  defect  are  dangling  bonds  (or, 
more  generally,  weakly  bonded  states'.  The 
number  of  these  states  is  considerably  smaller  than 
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expected  because  the  random  lattice  undergoes  re¬ 
laxation  processes  around  a  vacancy  (or  a  cluster  of 
vacancies)  which  result  in  a  substantial  bond  recon¬ 
struction  and,  consequently,  partial  elimination  of 
dangling-bond  states.  This  “healing”  process  can¬ 
not  be  complete  because  of  the  high  coordina¬ 
tion13, 19  of  Si  which  imposes  severe  geometrical  re¬ 
strictions  on  the  reconstruction  process. 

Calculations  for  the  ideal  single-vacancy  prob¬ 
lem20-22  show  that  dangling  bonds  are  associated 
with  local  states  whose  eigenenergy  lies  in  the  gap 
(about  0.75  eV  above  the  valence-band  edge).  Al¬ 
lowing  the  atoms  around  the  ideal  vacancy  to  relax 
pushes  these  states  toward  the  valence  and  the  con¬ 
duction  band.  On  the  basis  of  these  calculations, 
one  can  conclude  that  the  dangling  bonds,  which 
have  survived  reconstruction,  create  states  in  the 
gap.  These  states  are  of  the  order  of  1019— 1020 
cm-3eV-1  as  electron-spin-resonance  experiments 
show.11,12  It  is  generally  believed  that  hydrogen 
passivates  the  dangling-bond  states  by  forming 
Si  — H  bonds  which  are  associated  with  states  lying 
well  within  the  valence  band.  Furthermore,  if  the 
material  is  grown  in  the  presence  of  hydrogen, 
much  of  the  reconstruction  is  prevented  as  a  result 
of  Si  —  H  bond  formation.3  This  explains  why  the 
highest  amount  of  hydrogen  in  a-SiHx  is  so  much 
larger  than  the  number  of  dangling  bonds  (which 
survived  reconstruction)  in  the  unhydrogenated 
specimens.  One  may  view  the  abundant  presence 
of  hydrogen  during  the  growth  process  as  effective¬ 
ly  reducing  the  coordination  of  the  resulting  struc¬ 
ture  and  thus  allowing  the  growth  of  an  un¬ 
strained,  chemically  stable  substance.  The  widen¬ 
ing4-3  of  the  gap  upon  hydrogenation  has  been  at¬ 
tributed  to  the  stronger  Si  —  H  bond  as  compared 
with  the  Si  — Si  bond.  Here,  as  well  as  in  a  prelim¬ 
inary  report23,24  of  this  work,  we  argue  that  an 
equally  important  contribution  comes  from  the  ef¬ 
fective  reduction  of  the  pptr  interaction  upon  hy¬ 
drogenation. 

Theoretical  models 

From  the  above  discussion  it  follows  that  a  com¬ 
plete  theory  of  hydrogenated  a-Si  has  to  deal  with 
the  following  aspects  of  the  problem. 

(a)  Topological  disorder.  Models  incorporating 
topological  disorder  (TD)  using  a  continuous  ran¬ 
dom  network,  have  been  used25,24  in  conjunction 
with  the  orthogonalized-linear-combination-of- 
atomic-orbitals  method  to  perform  calculations  of 
the  electronic  structure  of  a-Si.H.  These  calcula¬ 


tions  are  consistent  with  photoemission  experi¬ 
ments,9  but  they  do  not  seem  capable  of  obtaining 
detailed  information  for  the  densities  of  states 
(DOS)  in  the  gap  region  and  in  the  conduction 
band.  Recently,  Cohen  et  at.11  argued  that  TD 
widens  the  gap  and  creates  a  tail  in  the  DOS 
which  enters  the  gap.  Using  small  Si-H  molecules 
terminated  by  Si  Bethe  lattices,  Allan  and  Joanno- 
poulos23  have  examined  the  question  of  ring  statis¬ 
tics  and  its  influence  on  certain  regions  of  the  spec¬ 
trum.  We  feel  that  increasing  hydrogenation  may 
reduce  the  importance  of  ring  statistics.  Finally, 
TD,  introduced  by  allowing  variation  in  the 
dihedral  angle,  seems  to  effectively  reduce  the  size 
of  the  ppir  interaction23  and  thus  contributes  fur¬ 
ther  to  the  widening  of  the  gap. 

(b)  Reconstruction.  Reconstruction  is  probably 
the  most  important  aspect  of  unhydrogenated  a-Si. 
A  qualitative  lattice  distortion  model  has  been  pro¬ 
posed  by  Watkins.29  This  model  has  been  used  as 
the  basis  to  address  the  problem  of  the  single 
reconstructed  vacancy  in  Si  using  elaborate 
Green’s-function  techniques.30,31  Also,  White  and 
Ngai32  have  discussed  reconstructing  states  at  the 
Si-SiO*  interface.  However,  it  seems  that  the 
amount  of  reconstruction,  and  therefore  its  impor¬ 
tance,  is  reduced  when  hydrogen  is  present  during 
the  growth  process.3  Thus,  depending  on  the 
abundance  of  hydrogen,  the  method  of  preparation, 
and  other  details  of  the  growth  process,  the  impor¬ 
tance  of  reconstruction  may  vary  from  a  dominant 
role  to  an  insignificant  detail. 

(c)  Chemical  nature  and  statistics  of  the  hydrogen 
incorporation.  By  chemical  nature  we  mean  wheth¬ 
er  the  hydrogen  is  always  bonded  to  one  of  the  four 
spi  hybrids  of  Si,  or  may  participate  in  other  bond¬ 
ing  configurations.  Even  if  hydrogen  is  only  bond¬ 
ed  to  Si  one  still  has  to  know  statistical  informa¬ 
tion  such  as  the  percentage  of  monohydrides  versus 
polyhydrides,  whether  or  not  there  are  some  hydro¬ 
gen  clustering  tendencies,  etc.  Obviously,  these 
questions  affect  the  electronic  structure  of  the  ma¬ 
terial. 

With  the  exception  of  the  continuous-random- 
network  work,23,24  most  attempts  in  the  literature 
to  study  the  above  three  aspects  are  based  on  con¬ 
sidering  small  Si-H  molecules.  These  Si-H  clusters 
are  either  isolated  or  terminated  by  hydrogen33  or 
by  Si  Bethe  lattice23  to  avoid  unphysical  boundary 
effects.  These  approaches  are  very  useful  in  re¬ 
vealing  certain  qualitative  and  even  semiquantita- 
tive  aspects  of  the  subject;  they  are  also  necessary 
in  some  cases  (e.g.,  other  hydrogen  bonding  config- 
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urations)  given  the  difficulty  of  the  problem. 
However,  these  approaches  cannot  be  considered  as 
a  substitute  for  calculations  dealing  with  macro¬ 
scopic  size  systems.  We  should  not  lose  sight  of 
the  fact  that  these  calculations  more  properly  ap¬ 
ply  to  the  study  of  Si-H  molecules  than  to  the  Si-H 
solid.  In  addition,  the  Bethe-lattice  calculations 
are  based  on  a  first-nearest-neighbor  interaction 
Hamiltonian  which,  as  has  been  shown,"  34  is 
inadequate  to  reproduce  the  correct  gap  or  a 
reasonable  conduction  band. 

The  work  we  report  here  is  one  of  the  first  at¬ 
tempts  to  deal  with  the  problem  at  the  actual  mac¬ 
roscopic  scale  which  involves  many  particular  local 
configurations.  We  found  it  necessary  to  omit  as¬ 
pects  related  to  topological  disorder  and  recon¬ 
struction,  and  concentrate  our  efforts  on  the  third 
aspect  which  is  the  role  of  hydrogen  in  the  elec¬ 
tronic  structure  of  a-SiH* .  We  think  that  for  fully 
hydrogenated  samples  with  high  hydrogen  content, 
hydrogenation  is  the  most  important  aspect  and 
that  TD  and  reconstruction  will  change  our  results 
only  quantitatively.  To  treat  the  role  of  hydrogen, 
we  have  used  a  particular  model  of  hydrogen  in¬ 
corporation  which  assumes  that  all  hydrogens  are 
bonded  to  Si.  Although  our  model  may  not  be 
quite  realistic,  we  think  that  it  incorporates  the  im¬ 
portant  features  (which  are  independent  of  particu¬ 
lar  models)  such  as  bonding  and  antibonding  states 
made  out  of  Si spl  and  hydrogen  orbitals,  a 
stronger  Si  — H  bond  (as  compared  with  Si  — Si 
bond)  which  affects  the  states  at  the  bottom  of  the 
conduction  band,  and  an  effective  reduction,  with 
hydrogenation,  of  the  ppir  interaction  which  is  very 
important  for  the  states  at  the  top  of  the  valence 
band. 

Thus  we  believe  that  the  main  conclusions  of  our 
work  have  a  much  wider  validity  than  the  particu¬ 
lar  model  from  which  they  were  derived.  Recently, 
a  calculation  complementary  to  ours  was  reported 
by  Pickett.13  Pickett  has  employed  the  self-con¬ 
sistent  pseudopotential  method  with  a  supercell 
configuration  to  study  the  electronic  states  of  the 
hydrogen-saturated  vacancy  in  Si.  His  approach 
differs  from  ours  in  that  he  is  using  a  periodic  ar¬ 
ray  of  atoms,  while  we  are  using  a  random  array. 
The  results  of  the  two  calculations,  however,  have 
the  same  qualitative  features.  In  the  Appendix  we 
utilize  his  conclusion  of  a  strong  Si  —  H  bond  and 
his  charge-density  contours  to  estimate  certain  H- 
Si  matrix  elements  of  our  tight-binding  Hamiltoni¬ 
an.  Finally,  we  refer  to  a  calculation  along  similar 
lines  to  ours  reported  by  Divincenzo  et  al ,  'h  This 


calculation  deals  with  a  model  defect  which  is  an 
isolated  monovacancy  in  an  otherwise  perfect  crys¬ 
tal. 

The  rest  of  the  paper  is  organized  as  follows. 
Section  II  describes  our  mode!  configuration  of  hy¬ 
drogenated  a-Si;  Sec.  Ill  gives  the  theory  of  the 
coherent-potential  approximation  as  applied  in  the 
present  work;  Sec.  IV  discusses  the  results  and 
compares  with  experiment,  and  the  Appendix  deals 
with  the  estimation  of  the  H-H  and  H-Si  matrix 
elements. 


II.  THE  PRESENT  MODEL 

Our  model  describes  hydrogenated  Si  by  con¬ 
structing  an  effective  lattice  whose  sites  may  have 
probability  c  of  being  vacant  and  probability  1-c  of 
having  a  Si  atom.  In  addition,  we  have  assumed 
that  hydrogen  atoms  may  be  located  along  the 
lines  connecting  a  vacant  site  with  its  nearest 
neighbors,  as  shown  in  Fig.  Kb).  Thus  we  have  in¬ 
cluded  in  our  model,  at  random.  Si  sites,  vacancy 
sites,  and  sites  that  have  one,  two,  three,  or  four 
hydrogen  atoms.  Using  this  model  of  disorder,  we 
have  used  a  tight-binding  form  of  the  coherent- 
potential  approximation37,38  (CPA)  to  perform  de¬ 
tailed  calculations  of  the  electronic  DOS. 

The  starting  point  of  the  present  calculations  is 
a  Slater-Koster  (SK)  Hamiltonian  H;  the  bases  are 
four  Si  orbitals  (one  |s)  and  three  | p))  which 
have  been  taken  as  orthonormal.  The  matrix  ele¬ 
ments  in  this  basis  have  been  chosen  in  such  a 
way  as  to  reproduce  the  band  structure  of  crystal¬ 
line  Si  (both  valence  and  conduction  band)  rather 
accurately."  Such  an  accuracy  is  necessary  in  or¬ 
der  to  study  dangling-bond  states.  As  we  have  dis¬ 
cussed  in  the  Introduction,  we  have  neglected  the 
TD  except  for  the  following  point.  Because  the 
dihedral  angle  (i.e.,  the  angle  which  determines  the 
orientation  of  the  three  bonds  which  are  attached 
to  one  end  of  a  given  bond  with  respect  to  the  oth¬ 
er  three  bonds  which  are  attached  to  the  other  end 
of  this  given  bond)  varies  in  a  disordered  structure 
in  a  range  which  starts  from  the  eclipsed  confi¬ 
guration  all  the  way  to  the  staggered  configuration, 
the  pptr  interaction  fluctuates.  As  we  have  men¬ 
tioned  before,  the  pptr  interaction  is  very  important 
because  it  controls  the  position  of  the  top  of  the 
valence  band.  In  the  periodic  case,  thepprr  in¬ 
teraction  is  equal  to  the  difference  l  =  Exy<  tttI 
-  Ex  x  (  7  y  t  )  =yt- )'?.  where  y„  =  <  6  H  ,  2  > 
and  =  (  6  ,  II  j  3  )  =  (  b  i  //  .4 )  and  !  2 ) ,  j  3  ) , 
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FIG.  1.  (a)  Two-dimensional  view  of  the  atom  confi¬ 
guration  for  Si  showing  the  sp }  orbitals,  (b)  Two- 
dimensional  view  of  the  atom  configuration  for  Si-H 
showing  the  replacement  of  one  Si  atom  by  four  H 
atoms. 


1 4),  and  |  6)  are  the  sp]  hybrids  shown  in  Fig. 
1(a).  In  the  random  structure  ( 6 1 H  \  3 
(b\H  1 4)  so  that,  in  general,  there  are  three  dif¬ 
ferent  matrix  elements  (6 1 H  |2),  (6 1  //  |3>,  and 
(6 1 H  1 4).  We  have  found  how  these  matrix  ele¬ 
ments  vary  with  the  dihedral  angle  and  we  have 
taken,  as  a  measure  of  the  ppir  interaction,  the 
average  of  the  maximum  of  the  three  differences  of 
one  of  them  from  the  mean  value  of  the  other  two. 
To  perform  the  average,  we  have  assumed  that  the 
dihedral  angle  has  a  uniform  (constant)  probability 
distribution.  The  effect  of  this  is  to  change  the 
matrix  elements  E^y  (777)  and  Exx  (777)  from 
their  crystalline  values22  1 .407  and  0.276  eV, 
respectively,  to  the  values  1.39  and  0.31  eV.  This 
reduces  the  ppir  interaction  A  from  the  value  1.131 
eV  to  the  value  1.08  eV.  This  leads  to  a  recession 
of  the  top  of  the  valence  band  and,  hence,  a  widen¬ 


ing  of  the  gap  from  Eg  =  1.0  eV  in  the  crystalline 
case  to  1.05  in  the  “amorphous”  state.  We  believe 
that  if  TD  were  taken  into  account  in  a  more 
rigorous  way,  it  would  lead  to  a  more  substantial 
widening  of  the  gap  and  to  a  tailing  of  the  states 
into  the  gap. 

Our  model  also  omits  reconstruction;  this  is  a 
very  serious  omission  in  the  cases  where  the  dan¬ 
gling  bonds  have  not  been  passivated  by  hydrogen. 
For  those  cases,  our  results  should  not  be  taken 
seriously  except  to  say  that  if  all  these  dangling 
bonds  were  present  there  would  be  substantial  DOS 
in  the  gap.  On  the  other  hand,  in  the  fully  hydro¬ 
genated  cases  we  expect  that  reconstruction  is 
minimal  and  thus  its  omission  in  our  model  is  not 
a  serious  shortcoming. 

The  way  hydrogen  is  incorporated  in  our  model 
is  shown  in  Fig.  Kb),  where  an  ideal  vacancy  has 
been  created,  resulting  in  four  sp 3  Si  dangling 
bonds.  Such  vacancies  can  accommodate  up  to 
four  hydrogen  atoms  as  shown  in  Fig.  1(b).  The 
result  of  replacing  a  Si  by  four  hydrogens  as  shown 
in  Fig.  Kb)  is  to  replace  the  Si-Si  matrix  elements 
[given  in  Table  I  of  Ref.  22  with  Exy  (777)  and 
Ex  x  (777)  changed  as  discussed  above]  by  H-H 
or  H-Si  matrix  elements.  In  the  Appendix  we  ex¬ 
plicitly  estimate  the  effective  matrix  elements  asso¬ 
ciated  with  the  configuration  of  Fig.  1(b).  They 
are  the  following : 

£,',(000)=  —8.72  eV,  £^(000)=  -  1.6  eV, 

£»<TTT>=*  eV. 

E^x ( 77 7 )  =  1 . 96  eV,  £;,(jjj)=0.<*eV, 

£;.x(ttt>=°-98  eV  • 

Actually,  the  second-  and  third-nearest-neighbor 
matrix  elements  will  be  affected  by  the  replacement 
of  a  Si  by  four  hydrogens.  We  assumed  that  this 
additional  modification  is  much  less  significant 
than  the  diagonal  and  first-nearest-neighbor 
changes  and  thus  we  have  omitted  it.  It  must  be 
pointed  out  that  the  £'(000)’s  given  above  are 
matrix  elements  between  fictitious  s  and  p  orbitals 
associated  with  the  four  hydrogens  shown  in  Fig. 
1(b)  (they  are  the  same  linear  combinations  of  the 
hydrogen  orbitals  j  D,  j  2' ),  |3'),and  14')  as 
the  actual  Si,  s,  and  p  are  of  the  corresponding 
four  hybrids;  see  Appendix)  and  the 
£'<777)’s  are  matrix  elements  between  the 
fictitious  four  hydrogen  s  and  p's  and  the  s  and  p's 
of  the  nearest  Si’s. 
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III.  COHERENT-POTENTIAL  APPROXIMATION 

As  we  have  mentioned  in  the  previous  section, 
we  have  used  the  tight-binding  (TB)  CPA  (Refs.  37 
and  38)  to  obtain  the  DOS.  Since,  in  addition  to 
the  configuration  shown  in  Fig.  Kb),  we  have  also 


considered  cases  where  one  or  more  of  the  hydro¬ 
gens  shown  are  missing,  we  had  to  generalize  the 
TB-CPA  to  handle  these  additional  configurations. 
The  CPA  condition  of  zero  scattering  on  the  aver¬ 
age  leads,  in  our  case  to  the  following  equation: 


([-c)ussi-GtuSi)-l+x0ul)(i-Gtu0r'  +  ^-  2  uu(\-G'Uur'  +  ^  2  Oun-Gtu,,rl 

*1*1  6  1=1 


+  ^7-  2  U3l(i-G,6f>/)_,+.t4{74(i-Gtt/4rl=o , 
4  .-1 


where  Gt  is  an  effective  Green’s  function  obtained 
from  the  corresponding  crystalline  Si  Green’s  func¬ 
tion  G  by  replacing  e,  and  ep  with  the  CPA  self¬ 
energies  2,  and  2p,  respectively.  From  Ref.  22  we 
have  that  e,  =£,,(000)=  —3.953  eV  and 
ep  =  EXX{000)  —  1.512  eV.  G(E)  is  a  4x4  diagonal 
matrix  with  matrix  elements  and 

where  CM  =  (Os  |  (£  —H)~l  |  Os  >  and 
Gn  =  (Ox  |  (£—//)“'  |  Ox  >,  where  jOs>and  |  Ox  ) 
are  the  s  and  px  Si  orbitals  at  the  site  0,  and  where 
U$„  the  Si  scattering  matrix,  is 

[et-2,  0  0  0  | 


ep 


value  for  silane  (Appendix).  Th^  matrix  — 2  is  a 
diagonal  matrix  like  (2)  but  without  the  e ,  and  ep. 
Un  is  a  matrix  corresponding  to  the  six  equivalent 
configurations  where  two  hydrogen  atoms  are 
present  with  probability  .x;/6, 


U,,=5 


7i  /2  0  0 

_  r’i  r’i  0  o  _ 


0  0  oo  0 
0  0  0  oo 


5-2  , 


and  similarly  for  U2 2,  t/;3,  f/;4,  U:i,  and  U:b. 
The  matrix  element  y2  =  —  1.78  eV  is  estimated  in 
the  Appendix.  U\<  is  a  matrix  corresponding  to 
the  four  equivalent  configurations  where  three  hy¬ 
drogen  atoms  are  present  with  probability  x ,  /4, 


Uv,  the  vacancy  scattering  matrix  is 

00  0  0  0 

_  0  00  •  0  0 

U“  =  0  0  00  0  •  (3) 

0  0  0  00 

U k  is  a  matrix  corresponding  to  the  four 
equivalent  configurations  where  only  one  hydrogen 
atom  is  present  with  probability  of  occurrence 
*  1  /4, 

y\  0  0  0 

—  _  0  00  0  0  _  _ 

"•>=5  0  0  00  0  (4) 

0  0  0  00 

and,  similarly,  for  U 12,  0 ,3,  and  U14.  5  is  the 
matrix  which  accomplished  the  orbital  trans¬ 
formation  (Appendix).  The  hydrogen  matrix 
element  y\  is  taken  equal  to  —3.38  eV  which  is  its 


71  72  72  0 

_  Yi  Y\  72  0  _  _ 
Un=S  .  ,  ,  .  5-2, 

72  7  2  7i  0 

0  0  0  oc 


and  similarly  for  C/32.  f’33,  and  C-;i.  L\  is  a  ma¬ 
trix  corresponding  to  the  case  of  four  hydrogen 
atoms  present  with  probability  x4: 

^  7i  72  7:  72 

_  _  72  7i  72  72  _  _ 

C/4=5  ,  .  ,  .  5-2.  (7) 

72  72  7i  7: 

72  72  7:  7i 

The  probability  of  occurrence  of  the  configuration 
shown  in  Fig.  1(a)  is  denoted  by  1  —  c:  the  probabil¬ 
ity  of  the  configuration  of  Fig.  Kb)  is  denoted  by 
.x4.  Obviously, 

C=X0+.X|  +X2  -KXj  -t-.x4  .  IS) 
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In  the  case  where  there  is  no  statistical  correlation 
among  hydrogens  the  probabilities  x,  (/  =0, ....  4) 
can  be  expressed  in  terms  of  the  quantity  x,  which 
is  the  ratio  of  the  total  number  of  hydrogens  over 
the  total  number  of  lattice  sites,  and  the  quantity  c: 


pressions: 

Gu^+U  —  efG|i  —  <,Gn  —  c(e,  —  e, ) 

— e,  +  e,eiG||  =0  , 
(12)* 

Gi2^-p  +( 1  —  —  ^p  ~  c^ep  ~  €p  ) 

—  €l>+€p€'pG22=0  , 


The  case  x  =  4c  represents  the  fully  hydrogenated 
case  where  x4  =c  and  x(  =0,  /  =0,  ....  3;  this  is 
the  case  where  the  reconstruction  effects  are  ex¬ 
pected  to  be  minimal  and,  consequently,  our  model 
to  be  more  realistic. 

Since  the  Green’s  function  G ,  is  a  function  of 
the  self-energy  2,  the  CPA  condition  [Eq.  (1)]  is 
very  complicated  to  solve,  even  numerically.  For 
this  reason  we  have  solved  Eq.  (1)  for  the  limiting 
cases  x  =0  (no  hydrogen)  and  x  =  4c  (all  vacancies 
saturated  by  hydrogen).  For  x  =0,  Eq.  (1)  com¬ 
bined  with  Eq.  (3)  reduces  to 

(1— c)C/Si(l  — G>£7Si)-1— cG>-l=0  .  GO) 

Utilizing  symmetry  Eq.  (10)  results  in  the  fol¬ 
lowing  two  scalar  equations: 

1,—e,  —  c/G\i(E,2,,2p) , 

(11) 

Zp=€p— c/Gj2(£,2„2<()  . 

Equations  (11)  are  solved  simultaneously  for  2, 
and  Ip  using  a  Newton-Raphson  iterative  pro¬ 
cedure.  The  Brillouin-zone  summation  necessary 
in  evaluating  GU,G22  was  done  for  240  k  points  in 
the  irreducible  zone. 

For  x  =4e,  Eq.  (1)  reduces  to  the  following  ex- 

_ I 

iVs.  =  li^£l(_imTr[(l-G,Usir,G,]!  , 

IT 


where  e,  =  £,'.,  (000)=  -8.72  eV  and 
ep=E‘x  x( 000)=  —  1.6  eV  are  the  effective  H-H  ma¬ 
trix  elements  (Appendix).  Equations  (12)  are  also 
solved  for  2,  and  Ip  by  iteration.  Since  both  sets 
of  Eqs.  (11)  and  (12)  are  complex,  we  found  it  com¬ 
putationally  efficient  to  separate  them  into  their 
real  and  imaginary  parts  and  do  the  computer  code 
in  real  arithmetic. 

For  the  intermediate  values  of  x,  we  have  em¬ 
ployed  an  average  r-matrix  approximation  (ATA) 
instead  of  the  CPA.  We  have  done  this  to  simplify 
the  computational  effort  since  the  additional  errors 
are  small,  and  because  the  cases  with  x  <4c  are  not 
very  realistic  due  to  the  reconstruction  that  takes 
place.  The  ATA- like  approximation  used  for 
0<x  <4c  is  the  following.  We  define  2,(£jc) 
i  =s,p  by  the  relation 

2 *(£;x)  =  1-^  2,(£;0)4- ■^-2,(£;4c)  , 

(13) 

where  2,(£;0)  and  2,(£;4c)  have  been  obtained 
from  the  CPA  described  above. 

Having  determined  2,  and  Ip  [using  the  CPA 
condition  (1)  for  the  cases  x  =0 ,4c  and  Eq.  (13)  for 
x  =5^0, 4c]  we  can  obtian  G,(E);  then  in  all  cases  the 
DOS  is  given  from  the  following  expressions: 


(14) 


—  2  |-ImTr[(l-G,U1,)-'G,]!  +  -1  £[ -taTrld-G^r'G,]] 

‘*1r  1  6ir  i  =  l 


+  -ZT  £  [-lmlT[(\-G,U2lr'G']\  +  ^-\ 
4ff  ,t|  IT 

where  the  notation  is  the  same  as  that  of  Eq.  (1). 

In  the  special  case  of  the  fully  saturated  vacancy, 
i.e.,  x  =4c,  x4=c,  and  Xj  =x2  =x j  =x„  =0,  the 
above  Eqs.  (14)  become 

iVSl  =  — —  Im Tr  (2-rH) — X—  G,  , 

*  fSi~?H 

(15) 


-ImTr[(l-GcU4)-|GJ|  , 


These  equations  can  be  further  reduced  to  give  the 
s-  and  p-like  components  of  the  DOS.  The  total 
DOS  is  of  course  the  sum  of  .Vs,  and  .VH.  The 
CPA,  as  described  above,  treats  the  diagonal  disor¬ 
der.  The  inclusion  of  off-diagonal  disorder  in  the 
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TB-CPA  is  computationally  very  complicated  and 
at  present  there  is  no  established  “best”  technique. 

In  our  calculations,  the  off-diagonal  disorder  was 
treated  within  the  virtual  crystal  approximation 
(VCA),  i.e.,  the  nearest-neighbor  matrix  elements 
were  replaced  by  their  averaged  values: 


c<.y'  i  2  i  ' 


Si-Si 

bij 


( 


)  + 


x  fSi-H,  l  l 

2  b'  i  1  2  2 


) 


i,j—s,x,y,2  .  (16) 

The  VCA  is  a  good  approximation  if  the  differ¬ 
ence  between  the  Si-Si  and  the  Si-H  matrix  ele¬ 
ments  is  small  in  comparison  to  their  average 
value.  This  condition  is  not  satisfied  for  all  matrix 
elements  and  so  there  is  a  need  here  for  improving 
our  present  calculational  techniques.  Given,  how¬ 
ever,  the  complexity  of  introducing  off-diagonal 
disorder  in  the  CPA  and  that  there  are  already 
other  uncertainties  in  our  model,  we  decided  in  the 
present  stage  to  work  with  the  VCA  for  the  off- 
diagonal  disorder.  We  have  used  the  simple  CPA 
i  no  cluster  extensions)  and  have  assumed  that  there 
is  no  statistical  correlation  among  the  various  con - 
figurations  discussed  before.  At  this  point  we  must 
note  that  it  has  been  proposed  that  vacancies  [see 
Fig.  1(b)]  tend  to  cluster  together  as  to  form  micro¬ 
voids  and  internal  surfaces.13  39  We  have  found  that 
this  clustering  effect  effectively  reduces  the  value  of 
c  and  tends  to  create  some  internal  surface  states 
which  make  a  small  contribution  to  the  total  DOS. 
Thus  the  vacancy  clustering  effects  can  easily  be 
incorporated  in  our  model  by  appropriately  reduc¬ 
ing  the  value  of  c.  Let  us  add  that  in  the  presence 
of  adequate  hydrogen  during  the  growth  process 
this  clustering  effect  may  not  occur. 


FIG.  2.  CPA  densities  of  states  in  the  gap  region  for 
the  Si  vacancy  case  for  different  vacancy  concentrauons. 


sion  of  the  effects  of  reconstruction,  we  certainly 
overestimate  the  number  of  these  states  and.  there¬ 
fore.  we  don’t  attempt  any  quantitative  comparison 
with  experiment. 

Figure  3  deals  with  the  restoration  of  the  band 
gap  upon  hydrogenation.  Figure  3  corresponds  to 
c  =0.05  and  0<x  <4c.  It  is  evident  that  by  in- 


IV.  RESULTS  AND  DISCUSSION 

In  Fig.  2  we  show  the  DOS  in  the  gap  region  for 
the  configuration  Si  and  vacancies  with  no  hydro¬ 
gen  introduced  yet.  W’e  note  the  appearance  of 
dangling-bond  states  in  the  gap.  The  density  of 
these  states  increases  with  vacancy  concentration  c. 
Also  the  gap  becomes  smaller  with  increasing  c  un¬ 
til  it  is  completely  filled.  It  is  also  interesting  to 
note  that  the  gap  states  have  as  their  center  of 
gravity  the  position  of  the  bound  state  '0.75  eV)  of 
the  ideal  single  vacancy. -0-—  It  should  be  stressed 
here  that  the  results  of  Fig.  2  are  useful  in  demon¬ 
strating  qualitatively  the  existence  of  dangling- 
bond  states  in  the  gap.  However,  due  to  the  omis- 
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FIG.  3.  Densities  of  states  of  SiH,  in  the  gap  region 
for  different  hydrogen  concentrations  x  c  =0.05  1 
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creasing  the  hydrogen  content  x,  the  states  in  the 
gap  are  reduced  until  they  are  completely  eliminat¬ 
ed  for  the  fully  hydrogen-saturated-vacancy  case 
(x  —0.20).  This  is,  of  course,  what  the  experimen¬ 
tal  situation  is  and  the  reason  why  a-Si:H  has 
semiconducting  properties  similar  to  those  of  crys¬ 
talline  Si. 

In  Fig.  4  we  show  the  DOS  of  Si-H  (with 
c  =0.05  and  x  =0.20)  for  the  whole  spectrum,  in¬ 
cluding  both  the  valence  and  conduction  bands. 

We  note  first  a  band  gap  Et  =  1.4  eV  that  is  wider 
by  0.4  eV  than  the  corresponding  Eg  =  1.0  eV 
which  our  tight-binding  Hamiltonian  gives  for  the 
nonhydrogenated  case.  This  widening  of  the  gap  is 
due  to  a  narrowing  of  the  third  DOS  peak  at  the 
top  of  the  valence  band.  We  have  identified  this  to 
be  the  result  of  a  decrease  of  the  difference: 
Ex.y(,j-fj)—Ex.x{---),  which  is  a  measure  of 
the  ppir  interaction.  This  recession  of  the  valence 
band  by  0.4  eV  is  in  excellent  agreement  with  the 
photoemission  experiments  of  von  Roedem  et  al.9 
Figure  4  also  shows  the  DOS  decomposed  by  site. 
The  hydrogen-site  DOS  shows  pronounced  peaks 
at  5.2,  7.6,  and  13.5  eV  below  the  Fermi  level. 
Comparison  with  photoemission  data  is  more  ap¬ 
propriately  done  after  smoothing  the  H-site  DOS 


FIG.  4.  Total  and  site-decomposed  densities  of  states 
for  SiH,  with  x  =0.2.  Note  that  the  Si  and  H  DOS  are 
multipled  by  ( 1  — c)  and  c,  respectively  (c  =0.05).  The 
Fermi  level  is  located  in  the  middle  of  the  gap. 


ENERGY  lev) 

FIG.  5.  Broadened  H-site  dpnsity  of  states.  The  ar¬ 
rows  indicate  the  peak  positions  from  the  photoemission 
measurements  (Ref.  9). 


by  applying  a  Lorentzian  broadening.  This  is  plot¬ 
ted  in  Fig.  5  which  shows  that  the  5.2  and  7.6  eV 
peaks  are  predicted  in  almost  exactly  the  same  po¬ 
sition  found  in  the  measurements.9  The  13.5  eV 
peak  is  not  seen  experimentally  for  reasons  we  do 
not  understand. 

The  widening  of  the  band  gap  has  been  demon¬ 
strated  experimentally  by  a  variety  of  experimental 
techniques.6-*  To  compare  with  these  experi¬ 
ments,  we  have  performed  a  calculation  of  the  joint 
DOS  NJ[E).  The  calculated  [iV‘,(£)]l/2  is  assumed 
to  be  proportional  to  the  measured  quantity 
(at£)l/\  where  a  is  the  absorption  coefficient.  A 
comparison  with  the  measurements  of  Cody  et  al. 1 
is  shown  in  Fig.  6,  where  [^Vy(£)]l/:  has  been  nor¬ 
malized  to  the  value  of  (a£)1/J  at  £=4  eV.  The 
experimental  graph  was  obtained  at  approximately 


FIG.  6.  Plot  of  the  square  root  of  the  joint  density  of 
states  versus  energy.  The  line  joining  the  open  circles 
represents  the  measured  quantity  (a£),/2  (Ref.  7). 
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16%  hydrogen  content,  while  the  theoretical  graph 
corresponds  to  20%.  Despite  this  and  the  fact  that 
the  theory  predicts  a  smaller  gap  the  overall  agree¬ 
ment  is  rather  good. 

Our  calculations  can  also  be  used  to  compare 
with  Auger  and  soft-x-ray-emission  measure¬ 
ments.40  We  present  in  Fig.  7  a  decomposition  of 
the  Si-site  DOS  of  SiH  in  its  s  and  p  components. 
The  r-like  and  p-like  DOS  are  proportional  to  the 
K  and  L  spectra,  respectively. 

We  now  turn  to  a  discussion  of  hydrogen- 
induced  antibonding  states  at  the  bottom  of  the 
conduction  band.  Let  us  first  examine  the  case  of 
a  single  configuration  of  the  kind  shown  in  Fig. 

Kb)  embedded  in  a  Si  lattice. 

The  bound  states  around  this  four-hydrogen 
cluster  will  be  given  by 


G„(£)= 

1 

(17a) 

£;(000)-£„(000)  ’ 

G22<£)  = 

1 

(17b) 

£^(000)-£„(000) 

The  graphical  solution  of  Eq.  (17a)  is  shown  in 
Fig.  8.  As  we  see  from  Fig.  8,  there  is  no  real 


FIG.  7.  Si-site  density  of  states  and  its  s-  and  p-like 
components. 


solution  of  Eq.  (17a)  because  the  intersection  of 
ReG  with  \/E'  —  £  occurs  within  the  band  where 
ImG^O.  On  the  other  hand,  this  intersection 
occurs  near  the  bottom  of  the  conduction  band 
(CB).  The  physical  meaning  of  no  real  intersection 
is  that  there  are  no  true  bound  states  associated 
with  the  four-hydrogen  cluster  embedded  in  a  Si 
lattice.  The  fact  that  the  intersection  occurs  near 
the  bottom  of  the  CB  where  ImG  is  very  small 
means,  physically,  that  the  four-hydrogen  cluster 
creates  resonance  s  states  i.e.,  states  where  the  wave 
function  has  a  peak  around  the  cluster.  Thes* 
resonance  states  can  be  viewed  as  a  hybridization 
of  the  Si-H  s-antibonding  states  with  the  regular  Si 
states  at  the  bottom  of  the  conduction  band.  Such 
resonance  states  are  associated  with  a  lower  than 
the  regular  CB  mobility  (because  the  electron  is 
almost  trapped  around  the  hydrogen).  Evaluation 
of  Eq.  (17b)  showed  neither  bound  states  of  p 
character  nor  any  resonance  states  below  3  eV. 

The  suggestion  of  Moustakas  et  al., 10  bom  out 
of  their  photoconductivity  measurements  that  Si-H 
antibonding  states  form  at  the  bottom  of  the  con¬ 
duction  band,  is  strongly  supported  by  the  present 
calculations.  This  is  shown  in  Fig.  9  where  we 
have  plotted  the  ratio  1VH  /N,  as  a  function  of  E. 
Indeed,  this  graph,  in  addition  to  the  peaks  in  the 
valence  band  that  we  have  already  discussed,  shows 
a  pronounced  maximum  at  the  bottom  of  the  CB 
indicating  strong  H  participation  in  the  formation 
of  these  states. 

Finally,  we  will  discuss  the  variation  of  the  gap 
size  Et  with  hydrogen  content  x.  We  define  the 
Fermi  level  £f  =  (£c  +£„)/2  where  £„  is  the  top 
of  the  valence  band  and  Ec  the  bottom  of  the  con- 


FIG.  8.  The  real  part  of  the  s-Iike  Green's  function 
plotted  as  a  function  of  energy  for  the  case  of  the  single 
impurity  that  consists  of  four  hydrogen  atoms.  The  dot¬ 
ted  line  indicates  the  value  of  the  right-hand  side  of  Eq. 
(17a). 
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FIG.  9.  Ratio  of  H-site  density  of  states  .Yh  to  the 
total  density  of  states  .V,  for  SiH  tx  =0.2). 


duction  band.  The  variation  of  £„,  £c ,  and  £f 
with  x  is  shown  in  Fig.  10,  where  in  all  cases  we 
are  dealing  with  vacancies  hilly  saturated  by  hy¬ 
drogen,  i.e.,  x  =4c.  We  note  first  that  £e  is  essen¬ 
tially  constant.  This  is  due  to  a  cancellation  of  the 
effects  of  disorder  (tends  to  push  Ec  down)  by  the 
effect  of  a  stronger  Si  — H  bond  (tends  to  push  Ec 
up)  as  manifested  by  the  larger  value  of  the  param¬ 
eter  yj  (see  Appendix)  in  SiH,.  £„  depends  mainly 
on  the  matrix  element  difference  y6— y,  or 

111  I  I  I 

equivalently  Exy(  ~~~  )—£,.,(  777 )  which  is  a 
measure  of  the  pp-rr  interaction.  Hydrogen  de¬ 
creases  Yb—Ys  and  so  pushes  £„  down  and  widens 
the  gap.  As  a  result  EF  is  also  pushed  down. 

Our  work  shows  that  the  size  of  the  gap  depends 
essentially  on  two  parameters:  the  bond  strength  y3 
and  the  ppv  interaction  y6  — yj.  Hydrogen  incor¬ 
poration  effectively  increases  y3  and  decreases 
y6  — y,  thus  producing  a  wider  gap.  This  effect  of 


hydrogenation  on  y3  and  y6-y5>  being  essentially  a 
local  chemical  effect,  is  expected  to  transcend  the 
validity  of  our  present  model  and  thus  constitutes  a 
general  feature  of  hydrogen  incorporation  in  a  Si 
tetrahedral  structure. 

In  conclusion  many  important  properties  of  fully 
hydrogenated  a-Si  (such  as  widening  of  the  gap, 
Si-H  bonding  states,  Si-H  antibonding  resonance 
states)  depend  mainly  on  the  local  chemical  en¬ 
vironment.  Thus  these  properties  are  largely  in¬ 
dependent  of  the  particular  model,  provided  that  it 
satisfactorily  treats  the  limiting  case  of  nonhydro- 
genated  Si.  Our  present  model  satisfies  these  re¬ 
quirements  (with  the  exception  of  the  omission  of 
topological  disorder)  and  as  a  result  we  expect  our 
main  conclusions  to  have  a  much  wider  validity 
than  the  model  itself. 
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APPENDIX:  ESTIMATE  OF  HYDROGEN- 
HYDROGEN  AND  HYDROGEN-SILICON 
MATRIX  ELEMENTS 


FIG.  10.  Variation  of  the  bottom  of  the  conduction 
band  Ec  of  the  Fermi  level  EF  and  the  top  of  the 
valence  band  E„  versus  hydrogen  content  for  fully  hy¬ 
drogenated  samples. 


For  each  Si  we  have  the  four  orbitals  |s ),  |  x  ), 
| y),  | D  or  equivalently  the  foursp3  hybridized 
orbitals  shown  in  Fig.  1(a).  The  transformations 
from  the  one  set  to  the  other  are  the  following: 


ID 

ID 

1  2  > 

ID 

|3> 

=s 

ID 

|4) 

ID 

ID 

I'D 

ID 

|2> 

1  y) 

=s 

|3> 

ID 

|4> 
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for  the  Si  at  A.  For  the  Si  at  B  assuming  that  AB 
is  along  the  111  direction,  we  have 


15) 

IT) 

16)  _ 

1  X  > 

I  7)  ~5 

l/> 

|8> 

|T> 

IT) 

15) 

lx) 

16) 

— 

1  y) 

17) 

1?) 

18) 

where  ST  is  the  transposed  matrix  of  S.  The  4x4 
matrices  5  and  S  are  given  by 


1  1  1 
*4  i  -i 


i  i 
-i  -t 
i  -i 
-i  i 


_  i  i  -i  i 
s=\  i  i  -i 


i  -i  -i 
i  i  i 
i  -i  i 
i  i-i 


Following  the  SK  notation,  the  diagonal  and 
nearest-neighbor  matrix  elements  of  the  Hamiltoni¬ 
an  in  the  |  s  ),  |  x  >,  \y  ),  |  z  )  representation  are 
denoted  as  follows: 

£,.1(000)=<s  \H  |j>  ,  (A71 

£',  ,(000)  =  (je  \  H  \x) 

=  <y \H  \y)  =  (z \H \z)  ,  (A8) 

£».»<TTT>  =  <*  l#|T)  ,  (A9) 

« 

=  <?\H  |s>=  •••  ,  (A10) 

£*.<(777>=<*|tf  i*>  =  0rl"l>’> 

=  {l\H\z)  ,  (All) 

£:-.>-<TTT>=<^l  ^  l>'>  =  <^!^  U> 

=  <y\H \x)  =  {y\H  |z)  =  •  •  ■ 


Hirabayashi41  has  written  the  matrix  elements  of  H 
in  the  hybrid  representation  as  follows: 

n-(l|H|l>.  (A  13) 


2 

II 

* 

1  2)  . 

'A14) 

73=<5|  H 

n> . 

(A15) 

y*=<6|H 

ID. 

(A  16) 

Yt  =  (6\H 

13)  , 

(A17) 

Yt=(6\H 

1  2)  . 

(A  18) 

where  the  orbitals  1 6)  and  |  2}  are  along  an¬ 
tiparallel  directions. 

Using  Eqs.  (Al)  and  (A2)  one  easily  expresses 
the  y‘s  in  terms  of  the  fs  and  vice  versa.  We  have 

y,  =  <000)+3£„(000>],  (A19) 

J'2=7{£».i(000)-£,.,(000))  ,  (A20) 

T3=7[^777>-6£,.,(777) 

-3£„<yyy)-6Ex,,(}y7)]  ,  (A21) 

n=7l^.,<777’-2£».*(777> 

+  £x.x(yyy)  +  2£x.y(yyy)]  ,  (A22) 

*^[£m<777>  +  2£,.*<777> 

+£x,x( 777>~ 2£,.y( yy y>]  ,  (A23) 

r6=7[£„(777)  +  2£x,x(jjj) 

-3£*.x<7  7  7>  +  2£x.,( y-y-y  1]  ,  (A24) 

£„(000)  =  y,  +  3y2  ,  (A25) 

£x,x(000)=y,-y2  ,  (A26) 

^(777)=7(y3  +  6r4  +  6y5  +  3y6)  ,  (A27) 

^(777)=7(-y3-2y4+2y5  +  y6) .  (A28) 

£,x,*(777)=7(-y3  +  2y4  +  2y5-3yft)  ,  (A29) 

^./<777)=7(-r3  +  2y4-2y.  +  y6)  .  (A30) 

In  Fig.  1  (b)  the  Si  at  A  has  been  removed  and 
four  hydrogens  have  been  placed  as  shown,  in  or¬ 
der  to  passivate  the  Si  dangling  bonds.  Charge- 
density  contours  resulting  from  the  pseudopotential 
supercell  calculation  of  Pickett35  involving  the  con¬ 
figuration  shown  in  Fig.  Kb),  demonstrate  that 
there  are  hydrogen  p  and  even  d  components  in  the 
eigenfunction.  Furthermore,  these  charge-density 
contours  strongly  suggest  that  this  multiple  /  hy¬ 
drogen  state  can  be  approximated  by  an  s-only  or- 
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bital  which  is  displaced  towards  the  Si  site  by 
about  20%  as  shown  in  Fig.  1(b).  If  Op,  02-,  Oy, 
and  04-  are  the  positions  of  the  centers  of  the  dis¬ 
placed  hydrogen  s  orbitals  |  l’>,  1 2'),  |  3'),  and 
1 4'),  respectively,  the  distances  are  as  follows: 
BO,  =2.30  a.u.,  B02=5. 18  a.u.,  0|02=3.50  a.u. 

Comparing  Figs.  1(a)  and  1(b)  we  can  see  that 
the  role  of  the  four  sp3  hybrids  |  1 )  to  1 4)  is 
played  by  the  four  orbitals  |  1')  to  1 4').  The  re¬ 
sult  of  this  replacement  is  to  change  the 
r bV2 . 7*  to  . .  . .  y'6  given  by 


F s  to  the  new  values  E" s.  Actually,  the  second- 
and  third-nearest-neighbor  matrix  elements  should 
be  modified  as  well.  These  modifications  are 
difficult  to  estimate  and  are  not  expected  to  be  as 
important  as  the  changes  in  the  diagonal  and 
nearest-neighbor  matrix  elements.  For  these 
reasons  we  have  omitted  these  modifications. 

We  now  proceed  to  estimate  the  y”s.  The  quan¬ 
tity  /1  is  taken  equal  to  its  value  in  the  SiH4 
molecule41 

Y,  =  —  3.38  eV  .  (A40) 


yi=<l’|tf  1 1‘>  , 

(A31) 

rj-<r|tf|2'> , 

(A32) 

yi=(5\H  1  r> , 

(A33) 

in . 

(A34) 

74i2)=<5 1 H  \  2')  , 

(A35) 

yi  =  (6  i  /f  |  3')  , 

(A36) 

Yt  =  (6\H  |2*>  . 

(A37) 

A  complication  associated  with  the  configuration 
of  Fig.  1(b)  is  that  y'md  and  yi,2l  are  not  necessarily 
equal  as  in  the  configuration  of  Fig.  1(a).  Howev¬ 
er,  as  we  shall  see  below  the  difference  between 
yMU  and  y*2l  turns  out  to  be  small,  thus  we  can 
replace  these  matrix  elements  by  their  mean  value, 
i.e.. 


yn  1 1 


(A38) 


Before  we  proceed  with  the  estimation  of  the  values 
of  the  Y’ s  we  mention  that  we  can  introduce  ficti¬ 
tious  s  and p  orbitals  |s'>,  \x'),  \y')>  \z')  associ¬ 
ated  with  four  hydrogens  of  Fig.  Kb)  through  the 
relation 


|s'> 

1  r> 

\x’) 

I  2') 

\y') 

=S 

13') 

\z') 

|4'> 

.(A39) 


These  orbitals  allow  us  to  define  hydrogen 
associated  E" s  by  replacing  in  Eqs.  (A7)  —  (A  12) 
the  |s>,  | x>,  |y>,  \z)  by  |s'>,  |.v'>,  |/>,  |z'>. 
Because  the  transformation  in  Eq.  (A39)  is 
identical  to  that  in  Eq.  (A2),  it  follows  that  the 
E” s  are  given  in  terms  of  /’ s  as  in  Eqs.  (A25)  and 
(A30).  Thus  the  removal  of  a  Si  and  the 
placement  of  four  hydrogens  as  shown  in  Fig.  ltb) 
is  equivalent  to  changing  the  six  matrix  elements 


To  estimate  y2  we  need  to  obtain  the  off-diagonal 
matrix  element  between  the  orthogonalized  hydro- 
genic  wave  functions  associated  with  the  configura¬ 
tion  shown  in  Fig.  1(b).  Mattheiss43  has  examined 
this  problem  in  detail  for  a  system  of  six  hydrogens 
placed  in  the  comers  of  a  canonical  hexagon.  We 
think  that  the  nearest-neighbor  matrix  elements  do 
not  depend  so  sensitively  on  the  geometry  and  con¬ 
sequently  Mattheiss’s  results  can  be  used  10  obtain 
a  fair  estimate  of  y.  We  have  fitted  Mattheiss’s 
results  for  separations  R  =2,  3,  5  a.u.  with  a  qua¬ 
dratic  function  times  the  exponential  function 
[exp(  —/?)].  We  found  from  this  fitting  that 

Yz=  -27.07(1.491  -0.012R  +0.077/1 2  )e~n 

(A41) 

in  units  of  eV  which  for  R  =3.5  a.u.  gives 

y2=  — 1.78  eV  .  (A42) 

To  obtain  yj  we  shall  write  it  by  employing  Eq. 
(A3)  as  follows: 

yj  =  {(<.f|/f  |  l'>-3<jf \H  |  l'>>  .  (A43) 

Chadi44  has  found  that  the  nonorthogonalized  orbi¬ 
tal  | f,  )  to  which  | s)  reduces  as  the  overlap  goes 
goes  to  zero  is  proportional  to  /?exp(  —  1.04/?). 

The  fact  that  this  s  Si  decays  almost  exactly  like  a 
hydrogen  s  orbital  suggests  that  Eq.  (A41)  may  be 
used  to  obtain  (s  |  H  \  l’>.  However,  the  extra 
factor  of  R  in  |  )  would  cause  the  matrix 

element  to  decay  more  slowly  than  the  right-hand 
side  of  Eq.  (A41).  To  take  this  into  account  we 
write,  in  eV 

y 

(s  \H  ]  l’>  =  -  -rr^- 27.07(1.491  -0.072/? 
kh-h 

+0.077 Rz)e~R  , 


where 


(A44I 
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4>,.h  =  ~  |('/3/4+lV/3)(l+K) 

D  1 

+  -±^{2+R/4)  exp(  —R) 

is  the  off-diagonal  matrix  element  between 
exp(  —  R )  and  R  exp(  —  R )  and 

KH.H  =  -(1.5  +  1.5R+R2/6)exp(— fl) 

is  the  same  matrix  element  between  two  hydrogen- 
ic  wave  functions.  To  check  the  accuracy  of  Eq. 

( A44)  we  substitute  for  R  the  Si  —  H  distance  in 
SiH4(R  =2.8  a.u.)  and  we  obtain  —3.58  eV.  This 
is  in  surprisingly  good  agreement  with  the  esta¬ 
blished42  value  of  —3.57  eV.  To  obtain 
(x\H  |  1')  we  take  into  account  that  the 
nonorthogonalized  |  x„  )  is  proportional  to  x  )  s„  ). 
Hence  it  is  plausible  to  write 

(x\H  I  r>=-cRx<s  |//  I  r>  ,  <A45) 

where  Rx  is  the  .x  component  of  the  vector  0,1?, 
and  the  constant  c  can  be  determined  from  the 
known  values  of  <Jc  J  H  |  1')  and  (ft  H  \  l'>  for 
the  SiH4  molecule.42  We  found  c  =  y.  Substi¬ 
tuting  in  Eqs.  (A44)  and  (A45)  the  values 
R  =0,-Z?  =  2.3  a.u.  and  Rx  =  1.33  a.u.,  we  obtain 
<s \H  |  1')  =  —4.85  eV  and  <.x  | //  |  l')=2.15  eV. 
Thus, 

yj  =  -  5.65eV.  (A46) 

The  quantity  y'm)  can  be  written,  by  employing 
Eq.  (A3),  as 

nu>=7«T|//|i'>  +  <.x|//|  r» 

=  -1.35  eV.  (A47) 

By  employing  Eq.  (A3)  we  can  rewrite  the  rest  of 
the  y”s  as  follows: 

/«2)=T«T|tf  |2'>-<x  j/f  |2'> 

|2'>-(i-|/f|2'»,  (A48) 

yi={«s\H\y)-(x\H\y) 

+  (y\Hiy)  +  (l)H  !V»  ,  (A49) 

*/6  =  T(<i|tf|2'>-<x|ff  |2’> 

+  {y\H \T)  +  (z\H  |2'>) .  (A50) 

One  may  attempt  to  calculate  the  matrix  ele¬ 
ments  in  the  right-hand  side  of  Eqs.  (A48)  — (A50) 


by  using  Eqs.  (A44)  and  (A45).  However,  Eq. 

(A44)  will  overestimate  the  size  of  the  ( (s  |  H  !  2' ) ) 
matrix  element  because  of  the  presence  of  the  hy¬ 
drogen  at  1'  which  represents  an  effective  repulsive 
potential.  In  Mattheiss’s43  calculation  this  reduc¬ 
tion  of  the  next-nearest-neighbor  transfer  matrix 
element  is  so  large  that  the  matrix  element  be¬ 
comes  almost  zero  (actually  it  changes  sign  and  is 
positive).  In  the  present  case,  due  to  the  fact  that 
the  hydrogen  1'  is  closer  to  Si  and  that  there  is  the 
extra  R  factor  in  the  Si  wave  function,  one  expects 
a  smaller  reduction  of  about  50%  assuming  that 
only  half  the  space,  i.e.,  the  region  around  the  hy¬ 
drogen  at  2'  would  contribute  to  the  integral.  This 
crude  reasoning  suggests  that  '(s\H  j  2')  is  rough¬ 
ly  half  the  value  given  by  Eq.  (A44),  i.e., 

(s  \H  |  2’>  =  —0.42  eV  .  (A51) 

The  other  matrix  elements  entering  Eqs. 

(A48)  — (A50)  can  be  estimated  by  employing  Eq. 
(A45)  and  the  fact  that  the  vector  2'B  and  3 'B  are 
(1.327,  3.799,  3.799)  and  (3.799,  1.327,  3.799), 
respectively.  We  obtain  thus 

(yjH  \y)~(x\H  12'}  =0. 19  eV  ,  (A52> 

(x\H  |3’>  =  <z|//  |3')  =  0r \H \T) 


~{z\H  j  2')  =0.53  ev  ,  (A53) 

which  give 

*4j,=  _0.84eV  ,  (A 54) 

y'j  =  -0. 1 1  eV  ,  (A55) 

4=0.23  eV.  (A56) 

Using  Eq.  (A38)  we  obtain  for  y\ 

4  =  —  1. 10  eV  .  (A57) 


The  corresponding  E" s  are  obtained  from  Eqs. 
(A25)  — (A30), 

E;,s (000)  =-8.72  eV  , 

£;J,(000)=-1.60  eV  , 

E;.,<777>=-3-05  «V. 

£»..*  (  t  7  7 )  =  I  ■  96  eV, 

£;.x(777>=0.64eV, 

£;,l777>=0-98  eV- 
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The  self-consistent  pseudopotential  method  and  supercell  configuration  are  used  to  study  the  electronic  structure  of 
the  hydrogen-saturated  vacancy  |HSV|  in  Si.  A  unified  picf.re  is  obtained  by  comparing  three  systems:  |i)  bulk  Si  ino 
vacancy),  (ii)  the  vacancy  alone,  and  |iii)  the  HSV  complex  in  which  each  Si  dangling  bond  is  compensated  by  H  The 
vacancy  dangling-bond  stales  are  found  to  be  removed  from  the  gap  and  the  H-H  interactions  do  not  interfere  with 
the  formation  of  a  strong  H-Si  bond.  The  reappearance  of  the  gap  and  the  form  of  the  charge  density  in  the 
proximity  of  the  HSV  implies  a  stability  which  suggests  this  may  be  one  of  the  strain-reiieving  configurations  w  hich 
occurs  upon  hydrogenation  of  amorphous  silicon.  These  calculations  also  suggest  a  potentially  observable  accepior 
state  if  a  HSV-like  defect  can  be  created  in  crystalline  silicon 


L  INTRODUCTION 

The  recent  upsurge  in  interest  in- amorphous 
silicon  arises  primarily  from  two  sources.  The 
first  involves  the  possible  widespread  technologi¬ 
cal  application  of  these  materials,  particularly 
hydrogenated  amorphous  silicon  (<t-Si:H)  in  solar 
cells.  The  second  impetus  centers  on  answering 
fundamental  theoretical  questions  about  the  nature 
of  electronic  states  in  disordered  systems,  with 
the  existence  of  a  gap  and  the  structure  in  its 
proximity  being  of  central  importance. 

The  theoretical  description  of  disordered  sys¬ 
tems  is  complicated  greatly  by  a  lack  of  know¬ 
ledge  of  the  atomic  geometry.  The  problem  of 
disorder  on  an  atomic  scale  is  exacerbated  in  a 
a-Si:H  by  the  existence  of  voids  and  related  large- 
scale  inhomogeneities,1  i.e.,  disorder  on  a  mac¬ 
roscopic  scale  which  is  extremely  difficult  to  mod¬ 
el  theoretically.  Attempts  to  describe  the  effects 
of  the  lack  of  short-range  order  on  the  electronic 
spectrum  of  a-Si  and  a-Si:H  have  concentrated  on 
valence  force  field  (bond-stretching  and  -bending) 
effects  or  on  low-symmetry  H-Si  or  Si  dangling- 
bond  configuration.  Each  of  these  approaches  foc¬ 
uses  on  specific  features  of  the  actual  system  of 
interest  and  together  they  have  led  to  a  rudimen¬ 
tary  understanding  of  some  of  the  observed  phen¬ 
omena. 

One  potentially  important  area  which  has  seen 
very  little  theoretical  attention  is  the  careful 
description  of  H-Si  bonding  in  a  solid-state  en¬ 
vironment.  Obviously  low-symmetry  configura¬ 
tions  involving  H-Si  bonds  will  occur,  and  it  is 
particularly  in  such  low-symmetry  situations  that 
self-consistency  is  often  found  to  be  essential  in 
a  realistic  description  of  the  electronic  structure. 
The  present  work  was  motivated  by  these  consider¬ 
ations.  As  a  first  step  toward  an  understanding  of 
H-Si  bonding  in  the  solid  state,  the  (by  now  stand¬ 


ard)  self-consistent  pseudopotential  method  has 
been  applied  to  the  H-saturated  Si  vacancy.  Al¬ 
though  this  is  not  a  low-symmetry  configuration,  it 
provides  (as  will  be  shown)  a  basic  example  of  the 
removal  of  dangling-bond  gap  states  in  a  bulk  en¬ 
vironment,  much  like  what  has  been  found  to  oc¬ 
cur  at  surfaces  and  is  expected  by  hydrogenation 
of  a-Si.  The  method  that  is  used  in  this  study 
guarantees  that  the  bonding  electrons  react  to  a 
potential  which  is  itself  derived  from  the  elec¬ 
tronic  charge  density,  unlike  those  used  in  prev¬ 
ious  work  (discussed  in  Sec.  IV). 

In  Sec.  n  the  calculations!  details  are  reviewed 
briefly.  Section  n  is  devoted  to  an  exposition  of 
electronic  structure  of  the  H-saturated  Si  vacancy, 
and  comparison  to  both  the  idea!  vacancy  and  to 
crystalline  silicon  is  presented.  The  local  density 
of  states  is  used  extensively  to  provide  an  under¬ 
standing  of  the  "chemical  bonding”  of  H  and  Si  in 
an  extended  system.  Comparison  with  previous 
work  is  presented  in  Sec.  IV,  as  well  as  a  brief 
discussion  of  the  possible  connection  with  experi¬ 
mental  data.  A  summary  is  given  in  Sec.  V. 

II.  METHOD  OF  CALCULATION 
A.  Self-consistent  pseudopotential  method 

The  self-consistent  pseudopotential  (SCP)  method 
has  been  documented  elsewhere2  and  has  been  used 
in  numerous  applications.  In  the  present  study 
local  ionic  pseudopotentiais  which  have  been  util¬ 
ized  successfully  in  other  applications’-  *  were 
chosen.  The  pseudopotentials  were  fit  to  the  form 

v'iont<7)  =  ~ a[cos(67)  -  c'  exp(- -1q*) ,  (1) 

where  fl  denotes  the  atomic  volume  and  n(l  -  c ) 

=  l  to  yield  the  correct  17  — 0  limit.  The  constants 
b,  c,  and  -.1  for  Si  and  H  are  listed  in  Table  I. 

With  the  SCP  method  it  is  conventional  to  cal- 
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TABLE  I.  Ionic  pseudopotential  parameters  (in  a.u.) 
t>,  c,  d  in  Eq.  (1)  for  H  and  Si. 


Si4* 

H* 

b 

0.790  7 

0.2800 

c 

-0.352  0 

-1.5380 

d 

-0.01807 

-0.0070 

culate  the  Hamiltonian  matrix  elements  in  a  plane- 
wave  basis,  with  a  cutoff  plane-wave  energy  Ex. 

To  improve  convergence  of  the  eigenvalues,  ad¬ 
ditional  plane  waves  up  to  a  cutoff  E2  are  used  in 
the  Lowdin  perturbation  scheme.  In  this  study  £, 

=  3  Ry,  £,  =  6  Ry  were  found  to  give  convergence 
of  eigenvalues  to  about  0.1  eV,  which  is  sufficient 
for  the  present  purposes 

The  plane-wave  expansion  of  the  charge  density 
and  screening  potential  (for  which  the  Hedin-Lund- 
qvist5  local  density  functional  was  used)  was  car¬ 
ried  out  to  37  stars  (equal  to  485  plane  waves).  In 
the  early  stages  of  self-consistency  the  four  spe¬ 
cial  fc- point  set  was  used  in  calculating  the  charge 
density.  The  final  charge  density  was  self-con¬ 
sistent  on  the  regularly  spaced  ten-point  set. 

B.  Supercell  geometry 

A  cubic  cell  with  edge  length  a0=  10.2626  a.u. 
and  containing  eight  diamond- lattice  sites  is  used 
as  the  unit  cell  in  a  superlattice  configuration. 
Three  distinct  systems  have  been  studied.  As  a 
control  system,  the  band  structure  of  crystalline 
Si  was  calculated  with  the  use  of  the  8-atom  unit 
cell.  Next,  one  Si  atom  was  removed  and-the 
self-consistent  electronic  response  was  recalcu¬ 
lated,  without  allowing  atomic  relaxation.  This 
configuration  corresponds  to  a  regular  array  of 
12.5%  ideal  vacancies,  but  in  some  ways  (which 
we  discuss  in  Sec.  ni)  the  results  are  similar  to 
those  of  isolated  vacancy  calculations.  Although 
this  supercell  is, small,  interactions  between  va¬ 
cancies  are  transmitted  along  bonding  chains  for 
which  a  neighboring  vacancy  is  four  sites  away. 
Finally,  the  four  dangling  bonds  per  vacancy 
were  "compensated”  by  H  atoms  placed  symmet¬ 
rically  around  the  vacant  site.  The  H-Si  bond 
length  was  chosen  to  be  2.874  a.u.  (=  65%  of  the 
Si-Si  distance)  in  accordance  with  the  regular¬ 
ities  of  the  H-Si-Sij  system  discussed  by  Luco- 
vsky.” 

There  are  two  features  of  this  H-saturated  va¬ 
cancy  (HSV)  configuration  which  should  be  noted. 
First,  it  is  evident  that  there  will  be  interactions 
between  HSV’s  in  neighboring  cells.  This  limita¬ 
tion  is  less  serious  than  it  seems  initially,  since 
there  is  rather  little  interest  here  in  the  isolated 


HSV.  In  a-Si:H  there  will  be  nearby  H-Si  com¬ 
plexes  which  are  nearly  as  well  modeled  by  re¬ 
peated  HSV’s  as  by  any  other  specific  configura¬ 
tion.  (The  sharp  structure  in  the  density  of  states 
characteristic  of  crystalline  systems  of  course 
will  be  spurious.)  Secondly,  there  is  the  likeli¬ 
hood  of  H-H  interactions  at  a  given  HSV.  The  H-H 
distance  is  2.563  a.u.  and  the  H-H  interaction  is 
appreciable  but  hard  to  quantify,  although  the  lo¬ 
cal  density  of  states  which  we  calculate  provides 
some  indication  of-the  nature  of  the  H-H  inter¬ 
actions. 

C.  Local  densities  of  states 

The  wave  functions  and  associated  weights  (see 
below)  were  calculated  on  the  regular  35  fc-point 
grid  (vertices  of  64  tetrahedra)  in  the  irreducible 
Brillouin  zone.  The  total  and  local  densities  of 
states  (DOS)  were  calculated  with  the  tetrahedron 
method  of  Lehmann  and  Taut,1  in  which  both  en¬ 
ergies  and  weights  are  interpolated  linearly  within 
each  tetrahedron. 

For  (he  weights  we  have  used  the  fraction  of 
charge  of  the  state  which  lies  within  a  sphere  of 
radius  R  centered  at  various  sites  of  interest. 

The  chosen  sites  were  (1)  the  vacant  site,  (2)  the 
H  site,  (3)  the  H-Si  bond  center  (denoted  "H-Si 
bond”),  (4)  the  nearest-neighbor  Si  to  the  vacant 
site  f "Si(l )”],  (5)  the  Si-Si  bond  center  ("Si-Si 
bond”),  and  (6)  the  second-neighbor  Si  ["Si(2)”]. 

In  each  case  R  has  been  chosen  to  be  0.15a0  =  1.539 
a.u.,  which  is  slightly  more  than  half  the  H-Si 
bond  length.  Therefore  the  various  spheres  over¬ 
lap  considerably  and  the  sum  of  their  weights  is 
meaningless.  For  the  Si  vacancy  calculation  the 
same  notation  is  used  although  no  H  atom  is  ac¬ 
tually  present. 

III.  CALCULATIONAL  RESULTS 
A.  Charge  tiensi  ty 

A  contour  plot  of  the  self-consistent  pseudo- 
charge-densily  of  the  HSV  in  the  (110)  plane  is 
shown  in  Fig.  1.  Doth  Si-Si  and  H-Si  bonds  lie 
in  this  plane.  The  Si-Si  bonds  are  virtually  ident¬ 
ical  to  those  calculated  for  the  bulk  (not  shown) 
using  the  same  potential.  This  result  is  encoura¬ 
ging  as  it  indicates  that  interactions  between 
charge-density  disturbances  from  neighboring 
HSV's  are  rather  unimportant. 

The  peak  in  charge  corresponding  to  the  H-Si 
bond  is  also  evident  in  Fig.  1.  The  charge  density 
in  the  bonding  region  is  similar  to  that  calculated 
for  the  H-Si  bond  at  the  Sid  11)  surfaced'’  The 
large  value  of  charge  dpnsitv  at  the  vacant  Si  site, 

~  17  electrons  per  diamond-lattice  unit  cell,  sug¬ 
gests  significant  interactions  between  the  four  H 
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FIG.  1.  Contour  plot  in  the  (110)  plane  of  the  pseudo- 
charge-densit7  of  the  H-saturated  SI  vacancy.  The  units 
are  electrons  per  diamond-lattice  unit  cell.  Solid 
straight  lines  denote  SI— Si  bonding  directions  and  dashed 
lines  indicate  H-Si  bonds;  crosses  denote  vacant  Si 
sites.  For  clarity  the  abscissa  and  ordinate  lengths 
are  1.5  times  the  corresponding  supercell  dimension. 

atoms  surrounding  the  vacant  site.  However, 
apart  from  the  immediate  region  of  the  H  atoms 
the  charge  density  is  typical  of  Si,  particularly 
the  channel  of  nearly  zero  density  running  parallel 
to  the  bonding  chains. 

The  charge  density  near  the  Si  vacancy  (not 
shown)  was  found  to  be  very  similar  to  those  cal¬ 
culated  for  isolated  ideal  Si  vacancy.9  As  an  il¬ 
lustration  of  the  charge-density  deformation  which 
results  from  (conceptually)  placing  H  atoms  so 
as  to  terminate  the  vacancy  dangling  bonds,  the 
charge -density  difference  (HSV- vacancy)  is  shown 
in  Fig.  2.  Thi3  difference  in  density  is  altered 
from  that  arising  from  overlapping  spherical  H 
atom  densities  in  two  ways.  First,  there  is  the 
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FIG.  2.  Contour  plot  in  the  U10)  plane  of  the  chart;*. - 
density  difference  between  the  H-s  curated  vacancy  and 
the  ideal  Si  vacancy.  Notation  is  the  same  as  in  Fig.  1. 
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r  1G.  3.  Tuid  end  local  density  of  states  of  the  H- 
saturated  vacancy,  as  defined  in  the  text.  Note  that 
for  the  local  densities  of  states  the  large  peak  at  -  12.5 
eV  has  been  scaled  down  by  a  factor  of  10- 


afo  re  mentioned  interaction  between  H  atoms,  re¬ 
sulting  in  a  small  net  flow  of  charge  into  the  va¬ 
cant  site  (centroid  of  the  four  proton  sites).  Sec¬ 
ondly,  and  more  noticeable  in  Fig.  2,  there  is  a 
shifting  of  the  pseudocharge  maximum  from  the 
H  site  along  the  H-Si  bonding  direction.  This  ts 
primarily  a  result  of  the  dangling-bond  charge, 
which  is  known  in  the  vacancy9  to  "heal”  by  de¬ 
creasing  its  magnitude  in  the  (broken)  bond  posi¬ 
tion,  being  attracted  back  into  the  H-Si  bond. 

The  separaiion  of  the  charge-density  differences 
in  Fig.  2  again  is  a  testimony  to  the  small  inter¬ 
action  between  charge  disturbances  even  in  the 
small  supercell  being  used  here. 

U.  Density  of  stain 

The  total  density  of  states  (DOS)  at  the  HSV  is 
shown  at  the  bottom  of  Fig.  3.  The  overall  struc¬ 
ture  is  similar  to  that  of  crystalline  Si.  shown  at 
the  bottom  of  Fig.  4,  calculated  using  identical  k- 
point  sampling  and  Brilloutn-zone  integration  tech¬ 
niques  as  described  in  Sec.  II C.  Most  significant 
is  the  result  that  terminating  the  vacancy  dangling 
bonds  with  H  removes  from  the  gap  the  well  known 
ideal  vacancy  gap  states.9  which  in  the  present 
supercell  fill  in  the  gap  as  shown  in  Fig.  5.  This 
is  just  what  is  expected  in  a  chemical  bonding  pic¬ 
ture  if  the  H-H  interactions  within  an  HSV  do  not 
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FIG.  4.  Total  and  local  density  of  states  of  crystalline 
Si. 


cause  too  much  disruption.  The  calculated  gap 
for  the  superlattice  of  HSV’s  is  0.5  eV  compared 
to  the  crystalline  value  of  0.9  eV  (for  this  poten¬ 
tial  and  for  the  convergence  criteria  that  have 
been  imposed).  This  0.5-eV  “gap"  is  discussed 
more  fully  below.  The  only  qualitative  difference 
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FIG.  5.  Total  and  local  density  of  states  for  the  ideal 
Si  vacancy  (superlattice).  Cross  hatching  is  used  to 
emphasize  the  segments  of  the  partial  densities  of  states 
which  are  strongly  vacancy  related.  The  energy  zero 
denotes  the  top  of  the  filled  states. 


between  the  HSV  DOS  and  that  of  crystalline  Si  is 
the  splitting  off  of  a  band  in  the  lower- valence 
region  -13  «£  s  -  12  eV.  Comparison  to  the  DOS 
of  crystalline  Si  in  Fig.  4  suggests  that  in  the  case 
of  an  isolated  HSV  this  state  will  occur  as  a  reso¬ 
nance  at  the  bottom  of  the  valence  bands. 

Also  shown  in  Fig.  3  are  the  local  DOS’s 
(LDOS's)  as  described  in  Sec.  n  C.  For  compari¬ 
son,  the  corresponding  LDOS’s  for  the  Si  site  and 
Si-Si  bond  in  crystalline  Si  are  pictured  in  Fig.  4. 
(The  primary  difference  between  these  latter  two 
very  similar  LDOS’s  is  a  larger  density  of  bond¬ 
ing  states  within  4  eV  of  the  gap  in  the  Si-Si  bond 
LDOS.)  Figure  3  shows  that  the  split-off  bound 
state  at  ~-  12.5  eV  is  confined  to  the  region  of  the 
four  H  atoms,  i.e.,  the  vacant  site  and  the  H  posi¬ 
tions.  (Note  that  this  huge  LDOS  peak  has  been 
divided  by  ten  in  Fig.  3.)  The  other  energy  region 
of  interest,  that  is,  where  there  is  large  varia¬ 
tion  in  the  LDOS  within  the  snpercell,  is  centered 
on  the  fundamental  gap  0. 0-0.5  eV.  The  LDOS 
immediately  below  the  gap  is  primarily  due  to  the 
Si-Si  bonding  regions,  as  in  the  crystal.  The 
H-Si  bond,  H  site,  and  vacant  Si  site  show  a  pro¬ 
gressively  lower  density  of  bonding  states  in  the 
region  -  2  to  0  eV  which  indicates  the  electrons 
in  the  H-Si  bond  are  bound  more  tightly  than  in 
the  Si-Si  bond.  The  states  above  the  gap,  on  the 
other  hand,  are  strongly  H- related  states,  with 
the  lowest  unoccupied  state  at  T  being  localized 
(apparently)  and  centered  at  the  vacant  site.  This 
suggests,  for  the  isolated  HSV,  a  bound  state  near 
the  top  of  the  gap,  which  we  discuss  further  below. 
In  Fig.  5  the  total  and  local  DOS’s  for  the  vacancy 
superlattice  are  shown.  The  regions  of  states 
which  are  strongly  affected  by  the  vacancy  are 
crosshatched  for  emphasis.  No  gap  exists,  the 
bound  states  in  the  gap  of  the  isolated  ideal  Si 
vacancy  are  broadened  by  intervacancy  interac¬ 
tions  into  bands  more  than  1  eV  wide.  (We  note, 
however,  that  the  charge-density  disturbance  due 
to  the  vacancy  is  well  localized  even  for  this  small 
supercell.)  These  dangling-bond  states  lie  pri¬ 
marily  in  the  region  -0.5  to  5  eV,  although  some 
extend  down  to  -  2  eV.  In  addition,  the  s  state  on 
the  Si  nearest  the  vacancy  is  repelled  upward  into 
the  low  DOS  region  at  -  9  eV  by  the  effectively  re¬ 
pulsive  potential  of  the  vacancy.  These  results 
correspond  closely  lo  the  isolated  vacancy  results 
of  Ref.  9,  with  the  cross  hatching  at  -9  eV  in  Fig. 

5  corresponding  to  A,  symmetry  states,  and  cross 
hatching  at  0  to  1  eV  corresponding  to  the  Tt  sym¬ 
metry  gap  states,  in  this  reference. 

C.  II  related  stales 

The  similarity  of  the  energy  bands  of  a  super¬ 
lattice  of  HSV’s  (“Si-H,”)  to  those  of  Si  in  the  same 
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FIG.  G.  Energy  bands  along  the  <  1 1 1  >  lirection  of 
the  ouperlattice  of  IISV'3  i"SitH4”)  .md  o;  SuU  silicon 
for  the  same  unit  cell  (“SI,1' ).  Open  circles  at  T  in¬ 
dicate  the  hydrogen-related  states  which  are  discussed 
in  the  text  and  pictured  in  Fig.  7. 


unit  cell  (“Si,”)  is  evident  in  Fig.  6,  where  the 
bands  are  shown  with  the  valence -band  maxima 
aligned.  From  this  figure  it  is  possible  to  identify 
the  bulk  Si  states  which  are  strongly  perturbed  by 
the  HSV  complex.  In  essence,  this  makes  use  of 
the  fact  that  the  tetrahedral  arrangement  of  four 
hydrogen  atoms  in  the  HSV  is  isoelectronic  with 
Si  and,  in  fact,  it  is  not  too  dissimilar  in  many  of 
its  effects  on  the  overall  spectrum.  For  this  com¬ 
parison  it  is  expedient  to  focus  on  states  at  T, 
rather  than  at  other  positions  in  the  Brillouin 
zone,  where  many  characteristics  of  Bloch  states 
are  sensitive  to  supercell  size  and  geometry.  The 
strongly  perturbed,  and  therefore  H-related, 
states  at  T  occur  both  below  and  above  the  lower 
gap  at  -  12  to  -  11  eV  in  Fig.  3,  and  at  the  bot¬ 
tom  of  the  conduction  band  at  0.5  eV.  Contour 
plots  in  the  (110)  plane  of  the  charge  densities  of 


FIG.  7.  Contour  plots  in  the  (110)  plane  of  the  three 
hydrogen-derived  states  at  the  zone  center  of  the  H- 
saturated  Si  vacancy  superlattice.  The  charge  density 
is  normalized  to  one  electron  per  superlattice  unit  cell. 

these  states  are  shown  in  Fig.  7.  The  two  lower 
states  [Figs.  7(b)  and  7(c)]  are  derived  from  a 
symmetric  combination  of  the  four  H  Is  orbitals, 
with  the  splitting  arising  from  interaction  with  the 
Si  s  states  which  lie  in  the  same  energy  range. 
Lower- symmetry  combinations  of  the  H  Is  or¬ 
bitals  bond  with  the  Si  sp3  orbitals  throughout  the 
range  -  8  to  0  eV. 

The  lowest  “conduction-band”  state  at  0.5  eV 
also  is  derived  from  a  symmetric  combination  of 
H  orbitals  (presumably  2s  and  2 p).  The  small 
charge  density  of  this  state  in  both  the  H-Sl  and 
Si-Si  bonding  regions  indicates  a  nonbonding  (ra¬ 
ther  than  antibonding)  tendency  toward  the  Si  sp1 
orbitals.  This  behavior  suggests  that  this  state 
corresponds  to  a  bound  state  in  the  gap  for  the 
isolated  HSV.  The  supercell  used  in  this  calcula¬ 
tion  is  far  too  small  to  get  any  reasonable  esti¬ 
mate  of  the  energy  of  such  a  state,  that  is,  wheth¬ 
er  it  is  a  deep  or  a  shallow  acceptor,  if  indeed  it 
lies  within  the  gap.  Either  type  could  be  consistent 
with  the  charge  density  in  Fig.  7(a):  (1)  the  state 
is  deep  and  very  localized,  in  which  case  the  in¬ 
terstitial  charge  density  would  be  due  to  super¬ 
cell-induced  mLxing  with  the  conduction-band  states 
(which  have  large  densities  in  the  interstices),  or 
(2)  the  state  is  shallow  and  weakly  localized,  in 
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which  case  the  conduction-band  mixing  and  inter¬ 
stitial  charge  density  is  an  integral  component  of 
the  state.  Qualitatively,  the  large  density  within 
the  tetrahedron  of  H  atoms  tends  to  support  the 
deep  acceptor  point  of  view. 

IV.  DISCUSSION 
A.  Relation  to  previous  work 

Using  an  extended  HUckel  method,  Choo  and  Tong10 
investigated  an  HSV  in  Ge.  Their  system  con¬ 
sisted  of  a  47-atom  crystalline  cluster  with  sur¬ 
face  dangling  bonds  compensated  with  H  atoms  to 
remove  them  from  the  gap  region.  When  the  cen¬ 
tral  Ge  was  replaced  with  four  H  atoms  a  hydro¬ 
gen-related  acceptor  state  (with  40%  of  the  charge 
on  the  four  H  atoms)  appeared  in  the  upper  part 
of  the  gap,  much  as  appears  to  be  the  case  in  the 
present  study. 

Recently  Economou  and  Papaconstantopoulos11 
have  studied  the  electronic  structure  of  a  random 
array  of  HSV’s  in  Si  using  a  tight-binding  co¬ 
herent-potential-approximation  method.  Aside 
from  the  different  method  of  calculation,  their 
work  differs  from  the  present  work  in  essentially 
two  ways.  First,  they  considered  a  random  array 
of  5  %  HSV’s  compared  to  the  periodic  array  of 
12.5%  HSV’s  studied  here.  Secondly,  the  H-H  in¬ 
teraction  within  a  given  HSV  was  taken  to  be  the 
same  as  in  the  SiH4  molecule.  Their  H  LDOS 
near  the  fundamental  gap  is  very  similar  to  the 
“H-site”  curve  of  Fig.  3,  with  a  strong  H  com¬ 
ponent  in  states  immediately  above  the  gap  and 
a  comparatively  small  H  component  in  bonding 
states  just  below  the  gap.  The  first  H  LDOS  peak 
below  the  gap  occurs  somewhat  lower  (-  6  to  5  eV) 
than  in  the  present  study  (-  4  to  -  2  eV),  and  the 
splitoff  (by  the  H-H  interaction)  state  at  the  bot¬ 
tom  of  the  valence  band  is  broadened,  but  in  other 
respects  the  H  and  Si  LDOS’s  are  similar. 

There  are  several  other  calculations13  in  the 
literature  which  include  the  basic  H-Si-Si,  con¬ 
figuration.  However,  the  similarity  of  these  sys¬ 
tems  to  the  HSV  is  insufficient  to  warrant  a  de¬ 
tailed  comparison  with  the  present  study. 

B.  Relationship  to  expert/rental  systems 

The  basic  unit  in  the  HSV  studied  here  is  the 
bonded  H-Si  pair,  which  is  expected  to  occur  in 


a-Si:H,  no  doubt  in  several  different  types  of 
complexes.  Moreover,  the  HSV  is  one  of  the  sim¬ 
plest  configurations  which  exhibits  hydrogen  clus¬ 
tering,  a  phenomenon  which  has  been  observed13 
in  a-Si:H.  Unfortunately,  the  presently  available 
experimental  results  provide  only  very  indirect 
information  on  atomic  configurations.  The  LDOS’s 
of  Fig.  3  bear  some  resemblance  to  photoemission 
data14  on  unannealed  a-Si:H,  which  show  H-related 
peaks  at  6  and  11  eV  binding  energy.  The  calcu¬ 
lated  peaks  at  7  and  12.5  eV,  which  will  shift  if  the 
HSV  is  allowed  to  relax,  may  be  related  to  the 
observed  states.  However,  it  is  naive  to  expect 
much  quantitative  correspondence  between  the 
present  study  and  any  of  the  variety  of  experi¬ 
mental  samples  of  a-Si:H.  In  particular,  the  cal¬ 
culated  gap  of  0.5  eV  and  the  precise  location  of 
peaks  are  sensitive  to  the  supercell  size  and, 
therefore,  the  calculated  gap  bears  no  clear  re¬ 
lationship  to  the  (much  wider)  experimental  gap 
which  plays  a  central  role  in  device  performance. 

Nevertheless,  the  strong  H-Si  bond  and  the  res¬ 
toration  of  the  gap  suggest  that  the  (relaxed)  HSV 
may  be  a  stable  local  configuration  in  spite  of  the 
H-H  interaction.  Thus  the  HSV-type  complex 
could  beoneof  many  strain-relieving  configurations 
to  occur  upon  hydrogenation  of  a-Si. 

Finally,  it  should  be  noted  that,  if  indeed  an 
HSV- related  acceptor  state  does  occur  in  the 
crystal  (which  should  be  resolvable  on  a  strictly 
theoretical  basis  with  more  precise  calculations’), 
it  should  be  detectable  in  defected15  Si  or  Ge. 

Note  added  in  proof.  A  recent  self-consistent 
Green’s-function  study  of  the  isolated  HSV  (D.P. 
DiVincenzo,  J.  Bernhole,  M.  H.  Brodsky,  N.  O. 
Lipari,  and  S.  T.  Pantelides,  in  Proceedings  of  the 
Topical  Conference  on  Tetrahedrally  Bonded 
Amorphous  Semiconductors,  Carefree,  Arizona, 
1981,  in  press)  found  no  states  in  the  silicon  ener¬ 
gy  gap- 
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With  the  use  of  the  self-consistent  pseudopotential  method,  the  Heilman-Feynman 
theorem  is  applied  to  study  the  symmetric  relaxation  in  the  hydrogen-saturated  vacancy 
(HSV)  which  has  been  used  previously  to  model  the  electronic  properties  of  amorphous  sil¬ 
icon  hydride.  The  hydrogen  and  nearest-neighbor  silicon  atoms  are  found  to  relax  outward 
by  0.46  and  0.35  a.u.,  respectively.  The  primary  effect  of  this  relaxation,  which  is  driven  by 
the  Urge  H-H  interactions  in  the  ideal  HSV,  is  to  restore  the  depleted  local  H  density  of 
states  just  below  the  band  gap.  It  is  suggested  that  geometries  other  than  those  considered 


here  may  lead  to  more  stable  configurations. 

I.  INTRODUCTION 

Recently  several  investigators  have  used  the 
hydrogen-saturated  vacancy  (HSV)  in  Si  to  model 
the  electronic  structure  and  transport  properties  of 
hydrogenated  amorphous  silicon  (a-SiHx).  A 
tight-binding  coherent  potential  approximation 
(CPA)  treatment1  on  a  random  5%  “alloy”  of 
HSV’s  in  a  silicon  lattice  displayed  a  valence-band 
spectrum  very  similar  to  that  of  a-SiH*  as  well  as 
reproducing  the  experimentally  observed  widening 
of  the  gap  when  the  H  content  is  increased.  The 
self-consistent  electronic  structure  of  the  HSV  itself 
has  been  studied  from  both  the  supercell2  and  the 
isolated-defect3  approaches.  Both  of  these  studies 
found  the  density  of  states  (DOS)  immediately 
below  the  gap  to  be  strongly  depleted  around  the  H 
sites,  and  the  former  study  (Ref.  2,  hereafter  re¬ 
ferred  to  as  I)  also  found  a  DOS  increase  around  the 
H  sites  for  conduction  states  immediately  above  the 
gap. 

The  depletion  of  conduction-band  states  near  the 
gap  by  H  has  been  suggested3,4  as  the  cause  of  the 
low  hole  mobility  in  a-SiH*,  and  dc  transport  calcu¬ 
lations5  on  the  CPA  model1  have  borne  out  this 
qualitative  picture.  More  recently,  a  calculation*  of 
the  absorption  coefficient  (i.e.,  ac  conductivity) 
based  on  the  CPA  model  has  suggested  an  interpre¬ 
tation  of  the  difference  between  the  “optical  gap” 
and  the  “DOS  gap"  in  a-SiH,.  Thus  it  seems  that 
the  HSV  has  provided  a  reliable  basis  for  under¬ 
standing  several  of  the  crucial  properties  of  a-SiH*. 

The  interpretation  discussed  above  has  all  been 
based  on  the  electronic  structure  of  the  ideal  HSV. 
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However,  in  I  it  was  emphasized  that  the  large 
repulsion  between  H  atoms  in  the  ideal  HSV  was 
likely  to  lead  to  substantia]  atomic  relaxation  in  this 
defect,  and  a  preliminary  study  of  the  magnitude  of 
the  relaxation  has  been  reported.7  In  this  paper  we 
allow  the  symmetric  breathing  relaxation  of  the 
HSV  and  find  substantial  rearrangement  of  the  H 
local  DOS  both  above  and  below  the  gap  which 
must  be  reconciled  with  the  data  on  a-SiHx  if  the 
HSV  is  to  remain  a  viable  model  for  this  system. 

The  plan  of  the  paper  is  as  follows.  In  Sec.  II  the 
numerical  procedure  which  is  used  to  calculate  the 
forces  on  the  atoms  is  presented.  The  numerical  re¬ 
sults  and  their  interpretation  are  given  in  Sec.  III. 
Tht  last  section  is  devoted  to  a  discussion  of  the  im¬ 
plications  of  these  results  for  our  current  under¬ 
standing  of  the  electronic  properties  of  a-SiHx. 


II.  NUMERICAL  TECHNIQUES 

The  self-consistent  pseudopotential  (SCP)  method 
has  been  described  in  detail  elsewhere.'  The  Si 
and  H  local  pseudopotentials,  Hedin-Lundqvist 
exchange-correlation  potential,  supercell  size  corre¬ 
sponding  to  eight  Si  atoms,  and  calculation  of  local 
DOS  (LDOS)  are  as  described  in  I.  The  novel 
feature  of  the  present  work  is  the  calculation  of 
forces  on  atoms  using  the  Heilman-Feynman  (HF) 
theorem.8  Formal  aspects  of  the  HF  theorem  in 
density-functional  theory  have  been  discussed  else¬ 
where.9  However,  since  this  study  may  be  the  first 
of  its  kind  for  a  model  of  a  localized  defect  in  a 
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bulk  solid,  and  since  technical  aspects  of  the  calcu¬ 
lation  of  HF  forces  have  not  been  presented  in 
much  detail  elsewhere,  selected  details  of  the 
present  calculations  will  be  presented  in  this  section. 

A.  Heilman-Fey nman  force  calculations 

For  a  given  atomic  configuration  (  R; ) ,  the  force 
Fy  on  an  ion  at  R;  is  given  by 

Fj  =F^°  +  Fy  ,  (1) 

where  F jm  is  tie  direct  Coulomb  contribution  from 
all  other  ions,9  which  can  be  calculated  using._an 
Ewald  technique,10  and  F '  is  the  electronic  contri¬ 
bution  which  takes  the  classical  form 

Fj  =  f  d*r  p(r)Vuj(r  —  Rj)  (2) 

in  terms  of  the  local  ionic  pseudopotential  Uj.  The 
charge  density  p  should  be  the  “exact”  charge  densi¬ 
ty,  i.e.,  that  which  satisfies  the  local-density- 
functional  equations.  An  approximate  charge  den¬ 
sity,  which  might  lead  to  a  reasonable  total  energy 
Elol  j  R,  |  (which  is  variational  in  p)  and  thereby  to 
reasonable  forces  from  the  numerical  derivative 
A£,oj/ARy,  can  lead  to  widly  uncontrolled  results 
for  Fj  if  the  HF  expression  is  applied  directly.  An 
example  is  the  Gordon-Kim"  procedure  for  deriv¬ 
ing  the  forces  between  two  closed-shell  ions  from  a 
total  charge  density  approximated  by  the  sum  of 
two  spherical  ionic  densities.  Direct  application  of 
the  HF  theorem  to  this  approximate  charge  density 
leads  to  nearly  vanishing  forces  for  ionic  separa¬ 
tions  of  interest,  resulting  solely  from  the  exponen¬ 
tial  tail  of  one  ionic  charge  overlapping  the  nucleus 
of  the  other  ion,  whereas  the  Gordon-Kim  pro¬ 
cedure  of  calculating  Etot  in  an  intermediate  step  is 
known  to  be  quite  reliable.  _ 

Physically,  the  HF  force  F*  originates  from  the 
distortion  (dipolar  with  respect  to  R j)  of  the  valence 
charge  density  (away  from  that  of  overlapping 
spherical  atomic  densities)  due  to  bonding.  In  prac¬ 
tice,  the  calculated  (and  therefore  approximate) 
charge  density  must  be  ( I )  converged  with  respect  to 
k-point  sampling,  (2)  self-consistent,  and  (3) 
represented  sufficiently  generally  to  allow  the  im¬ 
portant  bonding-charge  distortions.  A  plane-wave 
representation,  either  of  the  total  (valence  pseudo-) 
charge  density  p  las  done  here)  or  of  the  distortions 
of  p  due  to  bonding,  seems  to  be  a  particularly  satis¬ 
factory  representation  as  it  is  not  fixed  to  particular 
atomic  positions. __ 

.For  the  HSV,  F '  was  calculated  from  the  expres¬ 


sion 

F'=«2G^(|G|)p(G)e",GR>,  (3) 

c<C 

where  VjiG)  and  pi G)  are  Fourier  coefficients  of 
and  p.  The  summation  limit  G  must  be  extended 
until  convergence  is  obtained.  Figure  1  shows  the 
convergence  of  F*  for  both  H  and  Si  as  G  is  in¬ 
creased.  The  important  contributions  to  F*  arise 
from  wavelengths  (2n-/G)>  a0/4,  where  a0 
=  10.263  a.u.  is  the  Si  lattice  constant.  For  all  re¬ 
sults  quoted  here  the  cutoff" was  taken  as  G2  =  21 
(2ff/a0)2=  10  Ry. 

Accurate  values  of  F^  require  a  charge  density 
which  is  more  accurate  than  is  necessary  for  simply 
determining  a  self-consistent  potential  or  the  total 
energy.  This  in  turn  requires  more  precise  eigenvec¬ 
tors  of  the  Hamiltonian  matrix  i/g  g..  A  common 
procedure7  in  the  SCP  method  for  obtaining  eigen¬ 
values  efficiently  is  to  construct  //gg.  matrix  ele¬ 
ments  up  to  a  large  cutoff  G\  =  E 2  and  use  the 
Lowdin  procedure12  to  fold  down  the  eigenvalue 
problem  to  a  smaller  one  corresponding  to  G]  =  Et. 
For  the  present  calculations  the  Lowdin  unfolding 
procedure12  for  subsequently  obtaining  the  eigenvec¬ 
tor  components  up  to  G=G2  has  also  been  em¬ 
ployed.  These  procedures,  which  are  reminiscent  of 
second-order  perturbation  theory  in  the  matrix  ele¬ 
ments  Wg  5  ■  for  Ei  <G2,  G'1  <E z,  produce  eigen¬ 
vectors  of  sufficient  accuracy  for  the  present  pur¬ 
poses.  The  cutoffs  £,=4.5  Ry,  E2  = 8  Ry  corre- 
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FIG.  1.  x  component  of  the  electronic  contribution  to 
the  force,  in  hartrees/bohr,  on  H  and  nearest-neighbor 
Si  3toms,  illustrating  convergence  vs  G-space  cutoff  <7. 
Results  are  for  configuration  1. 
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spond  to  plane-wave  basis  sets  of  approximately  175 
and  400  plane  waves,  respectively,  and  provide  only 
slightly  less  accurate  solutions  to  the  400  x  400  sec¬ 
ular  equation  at  substantially  less  cost  than  by 
direct  methods. 

In  Table  I  examples  of  the  changes  in  the  calcu¬ 
lated  forces  due  to  large-G  components  of  the  wave 
function  are  displayed.  For  this  convergence  check , 
the  single  special  k  point13  was  used  to  approximate 
the  charge  density  for  the  idea!  configuration,  and 
the  Lowdin  unfolding  procedure  was  applied  to  give 
initially  300  and  then  400  plane-wave  coefficients  of 
the  wave  functions.  The  calculated  forces  on  the  H 
and  Si  atoms  increased  in  magnitude  by  35%  and 
6%,  respectively,  while  the  total  energy  per  atom 
changed  only  from  —5.5028  to  —5.5058  Ry. 

The  set  of  four  special  k  points13  for  the  simple 
cubic  lattice  has  proven  adequate  for  the  HSV  su¬ 
percell  for  determining  p.  One  special  point  is 


26 

clearly  insufficient,  while  results  almost  identical  to 
those  from  the  four-special-point  mesh  were  ob¬ 
tained  using  the  10-  and  35-point  regular  meshes 
and  the  20-special-point  set. 

B.  Geometry 

The  HSV  supercell  is  a  cube  corresponding  to  the 
crystalline  Si  cubic  cell  of  edge  length  a0  and  con¬ 
taining  eight  Si  sites.  For  the  HSV  the  Si  atom  at 
the  origin  is  replaced  by  hydrogens  at  ( u,u,u)  and 
equivalent  positions  tetrahedrally  placed  around  the 
origin.  For  the  ideal  HSV,  u  was  chosen2  as 
0.0883a0.  corresponding  to  a  H  — Si  bond  length  of 
Sij-SiH  suggested  by  Lucovsky.14  Only  symmetric 
“breathing-mode”  relaxations,  i.e.,  atoms  relaxing 
radially  from  the  center  of  the  defect  only,  of  the 
*  HSV  have  been  considered.  Thus  only  u  and  the 


TABLE  I.  Comparison  of  the  x  component  of  the  calculated  electronic  forces  lin  arbitrary 
units)  on  H  and  Si  atoms  when  300  or  400  plane-wave  coefficients  are  included  in  the  plane- 
wave  expansion.  Contributions  from  each  star  of  reciprocal-lattice  vectors  as  well  as  the  total 
are  given.  Owing  to  the  symmetric  placement  of  the  Si  atom  {at  ( a,j/4,a0/4,a0/4)]  within  the 
unit  cell,  several  stars  G  give  vanishing  contributions  to  the  force  independently  of  piG). 


H  Si 


a 

2  TT 

G  ‘ 

300 

400 

300 

400 

100 

-1.57 

-1.60 

-9.54 

-9.70 

110 

'  -3.52 

-3.56 

111 

-8.78 

-8.78 

200 

-1.21 

— 1.23 

210 

-4.69 

-4.76 

1.64 

1.66 

211 

-1.49 

— 1.48 

-0.59 

-0.61 

220 

2.85 

3.13 

221 

1.61 

1.71 

1.71 

1.75 

310 

0.32 

0.31 

311 

11.86 

12.66 

222 

4.52 

4.68 

320 

-0.56 

-0.63 

2.68 

3.03 

321 

0.69 

0.77 

-4.68 

-5.26 

400 

2.52 

2.91 

322 

1.42 

1.58 

-1.54 

-1.73 

330 

0.39 

0.44 

-0.33 

-0.43 

331 

0.19 

0.18 

420 

-0.01 

-0.02 

421 

-0.01 

-0.01 

-0.05 

-0.07 

332 

0.06 

0.07 

0.61 

0.69 

422 

-0.04 

-0.05 

430 

0.21 

0.27 

-0.58 

-0.75 

431 

0.19 

0.24 

0.25 

0.38 

333 

0.06 

0.07 

Total 

5.01 

6.90 

-  10.42 

—  11.04 

-10.42 


-11.04 


26 


SYMMETRIC  RELAXATION  OF  THE  HYDROGEN-SATURATED  . . . 


5653 


coordinate  v  of  the  neighboring  Si  atom  at  (u,v,u) 
and  symmetric  positions  are  allowed  to  vary 
(u=0.25a0  f°r  the  ideal  HSV).  Changes  in  unit¬ 
cell  volume  have  not  been  considered. 


III.  RESULTS 
A.  Relaxation  of  the  HSV 

In  addition  to  the  ideal  HSV  geometry,  three  re¬ 
laxed  geometries  have  been  studied.  Their  coordi¬ 
nates  and  the  calculated  atomic  forces  are  presented 
in  Table  II.  Configuration  3  was  determined  from 
the  previous  configurations  (0,1,2)  by  fitting  forces 
determined  from  the  harmonic  potential-energy  ex¬ 
pression,  valid  for  small  enough  displacements, 

V(u,v)=Au  +Bv+Cu2+Dv2  +  Euir ,  (3) 

to  the  six  calculated  forces,  and  then  determining 
the  equilibrium  values  of  u  and  u.  After  calculating 
the  forces  for  configuration  3  (Table  II),  forces  de¬ 
rived  from  expression  (3)  were  least-squares-fitted 
twice,  to  the  calculated  forces  for  all  four  configu¬ 
rations,  and  separately  to  configurations  1,  2,  and  3 
(those  nearest  equilibrium).  From  these  two  poten¬ 
tial  functions  the  equilibrium  values  u0  and  v0  and 
their  uncertainties  (arising  from  which  fit  one 
chooses  to  use)  were  found  and  are  given  in  Table 
II. 

The  results  of  the  geometrical  relaxation  can  be 
summarized  as  follows.  The  H  atoms  relax  out¬ 
ward  by  1  u0  —  u(0)  |  =0.46  a.u.  [u(0)  =  u  for 
configuration  0]  while  the  neighboring  Si  atoms  re-  ' 
lax  outward  by  0.35  a.u.,  resulting  in  a  7%  decrease 
of  the  H  — Si  bond  length.  This  relaxation  rotates 
the  Si  — Si  bond  by  4.5*  and  compresses  it  by  2.5%. 
The  ideal  and  relaxed  geometries  in  the  (110) 
plane  are  shown  in  Fig.  2.  In  terms  of  u=u—u0. 


HSV  GEOMETRY 

INTERSTITIAL 

- IDEAL  SITE 


FIG.  2.  Schematic  geometry  of  the  HSV  relaxation  in 
the  (110)  plane  of  the  supercell.  Dashed  lines  represent 
the  configuration  before  relaxation. 


v—v  — 1>0.  the  potential  energy  V  for  small  displace¬ 
ments  can  be  written 

V(u,v)=au2  +  bv2  +  cu  v  .  (4) 

From  this  expression  the  energy  A E  gained  by  re¬ 
laxation  is  found  to  be  about  1  eV/supercell,  i.e., 
per  HSV  defect.  Direct  calculation  of  A£  (dis¬ 
cussed  more  fully  below)  between  configurations  0 
and  3  is  1.1  eV. 


B.  Changes  in  eigenstates  and  LDOS 

As  expected,  this  rather  large  atomic  relaxation 
leads  to  significant  changes  in  the  hydrogen-related 
states  and  LDOS.  In  Fig.  3  the  supercell  band 
structure  along  I~-R  is  shown  for  both  the  ideal  and 
relaxed  HSV.  The  states  of  particular  interest  in 
following  the  effect  of  relaxation  are  numbered  1 
through  6,  and  the  charges  of  each  of  these  states 
within  spheres  distributed  throughout  the  cell  are 
given  in  Tabic  III.  States  1—3  lie  in  the  lower 
valence  bands  and  before  relaxation  each  is  strongly 


TABLE  II.  x  components  of  H.u  and  Si,u  positions  (in  units  of  the  lattice  constant  10.263  a.u.),  and  the  x  com¬ 
ponents  of  ionic,  electronic,  and  total  forces  on  H  and  Si  atoms  (hartree/bohr).  Equilibrium  positions  are  calculated  by 
fitting  to  quadratic  potentials  as  described  in  the  text;  uncertainties  (in  parentheses)  are  due  to  the  two  separate  fits  to 
the  calculated  forces  and  do  not  include  uncertainties  in  the  calculated  forces. 


Configuration 

u 

V 

FT 

H 

F'< 

FT 

FT 

Si 

FI 

pi  r* 

*  t 

0 

0.0883 

0.2500 

-0.019 

0.038 

0.019 

0.060 

-0.04! 

0.019 

1 

0.0960 

0.2500 

-0.090 

0.091 

-0.001 

0.080 

-0.040 

0.040 

2 

0.1000 

0.2560 

-0.108 

0.115 

0.007 

0.036 

-0014 

0.022 

3 

0.1203 

0.2730 

-0.213 

0.202 

-0.01 1 

-0.049 

O 

o 

tjt 

u 

0.003 

Equilibrium 

0.1 1*2(3) 

0.2696(1) 

-0.179 

0.179' 

0.0 

-0.040 

0.040* 

0.0 

'Calculated  from  equilibrium  condition  FT=0. 
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localized  near  the  H  atoms.  The  outward  relaxation 
of  the  H  atoms  decreases  the  overlap  of  the  attrac¬ 
tive  H  pseudopotentials  in  this  region.  As  a  result 
these  states  rise  in  energy  and  lower  their  amplitude 
near  the  H  atoms,  although  state  2  (and  to  a  lesser 
extent  state  3)  remains  strongly  H  related. 

States  4— 6  lie  near  the  gap  and  therefore  are  cru¬ 
cial  in  determining  low-excitation-energy  properties. 
State  4  is  the  uppermost  valence  band  whose  eigen¬ 
value  is  used  to  fix  the  energy  zero  in  Figs.  3  and  4. 
[With  respect  to  the  average  potential  K,OI(G  =  0), 
the  energy  of  state  4  moved  downward  by  0.5  eV 
during  relaxation.]  From  Table  III  it  can  be  seen 
that  for  this  state  the  charge  near  H  nearly  doubles 
during  relaxation.  Conversely,  the  charge  on  the 
lower-conduction-band  state  5  decreases  dramatical¬ 
ly  near  and  between  the  H  atoms  during  relaxation, 
although  its  eigenvalue  changes  very  little.  Finally, 
.state  6,  the  low-conduction-band  state  at  R  in  the 
ideal  HSV,  is  raised  by  nearly  3  eV  by  the  relaxation 
although  its  strong  H-related  character  is  not 
changed  dramatically  (Table  III).  Both  states  5  and 
6  have  large  amplitudes  in  the  interstitial  regions  of 
the  cell. 

The  LDOS’s  for  the  spheres  along  the  bonding 
chain  pictured  in  Fig.  2  are  given  in  Fig.  4  along 
with  the  total  DOS.  The  only  notable  change  in  the 
total  DOS  occurs  at  the  bottom  of  the  valence 
bands,  where  the  strongly  H-related  peak  at  —  13  to 
- 12  eV  found  in  I  for  the  ideal  HSV  is  found  to 
merge  with  the  low  Si  valence  states  during  relaxa- 

TABLE  III.  Relative  charges  within  each  of  the  equal-volume  spheres  at  the  positions  indicated  in  Fig.  2  for  the  six 
states  designated  in  Fig.  3.  Both  ideal  (configuration  0,  in  parentheses)  and  near  equilibrium  (configuration  3)  values 
are  given,  with  normalization  such  that  a  uniform  state  will  have  a  value  1.00  within  each  sphere.  Also  given  below  is 
the  total  charge  (in  electrons)  within  each  sphere  and  the  ratio  "charge  (relaxed)Zcharge  (ideal)." 


State 

number 

Vacant 

site 

H 

H— Si  bond 

Sid) 

Si  — Si  bond 

Si(2> 

Interstitial 

1 

0.44 

1.13 

1.53 

1.61 

2.!5 

1.74 

0.61 

(7.90) 

(5.29) 

(2.72) 

(1.66) 

(1.06) 

(1.06) 

(0.14) 

2 

6.2-7 

7.92 

4.60 

0.73 

0.50 

0.60 

0.82 

(13.15) 

(7.35) 

(1.55) 

(0.54) 

(1.68) 

(1.15) 

(0.261 

3 

1.49 

3.06 

3.20 

2.45 

1.75 

0.38 

1.15 

(17.46) 

(11.00) 

(4.25) 

(1.29) 

(0.05) 

(0.25) 

(0.05) 

4 

0.66 

2.86 

2.20 

1.36 

2.75 

1.21 

0.05 

(0.83) 

(1.50) 

(1.52) 

(1.21) 

(1.40) 

(2.74) 

(0.07) 

5 

1.98 

1.10 

0.31 

0.63 

0.89 

1.37 

8.22 

(5.27) 

(2.18) 

(0.90) 

(0.77) 

(0.09) 

(0.54' 

(4  99) 

6 

8.64 

2.96 

1.61 

1.80 

1.12 

0.15 

4.59 

(8.58) 

(3.61) 

(2.13) 

(2.34) 

(0.18) 

(0.71) 

(1.571 

Total  charge 

0.84 

2.35 

2.21 

1.39 

1.34 

2.03 

0.32 

(2.10) 

(2.12) 

(1.90) 

(1.51) 

(1.79) 

!  1 .47) 

020) 

Ratio 

0.40 

1.11 

1.17 

0.92 

0.75 

1.38 

1.63 

- IDEAL. 

-  RELAXED 


FIG.  3.  Supercell  band  structure  along  the  T -R  direc¬ 
tion  for  configuration  0  (ideal)  and  configuration  3  (near 
equilibrium).  Dotted  arrows  indicate  states  which  are 
affected  strongly  by  the  relaxation.  Numbers  label 
states  discussed  in  text  and  referred  to  in  Table  III. 
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ENERGY  (eV) 

FIG.  4.  Total  and  local  densities  of  states  in  the  spher¬ 
ical  regions  pictured  in  Fig.  2.  In  the  LDOS  curves  the 
•dashed  peak  below  — 12  eV  has  been  divided  by  ten.  All 
LDOS  are  plotted  on  the  same  scale.  States  move  from 
the  single-  to  the  double-cross-hatched  regions  during  re¬ 
laxation. 


tion.  Accordingly,  the  vacant-site  and  H-site 
LDOS’s  show  a  general  decrease  in  the  — !  3-  to 
—  9-eV  region. 

There  is  a  large  increase  in  the  LDOS  just  below 
the  gap  at  the  H  site,  H  — Si  bond,  and  the  second- 
neighbor  Si  ["Si(2)"],  of -which  the  latter  may  be  a 
“supercell  effect”  which  indicates  that  for  these  en¬ 
ergies  the  charge  perturbation  extends  to  larger  dis¬ 
tances.  Conversely,  near  the  bottom  of  the  conduc¬ 
tion  bands  there  is  a  strong  decrease  in  the  LDOS  at 
the  vacant  site,  H  site,  and  H  — Si  bond,  together 
with  a  significant  increase  near  Si(2)  which  again 
may  be  a  supercell  effect. 

It  should  be  noted  that  little  can  be  learned  about 
the  energetics  of  the  relaxation  solely  from  the 
DOS,  i.e.,  from  the  sum  over  occupied  states  of  the 
eigenvalues.  It  is  found  that  the  change  in  eigen¬ 
value  [relative  to  TKjr,r„(G  =  0)=0]  sum  between 


configurations  3  and  0  is  —7.70  eV,  while  the 
change  in  the  Coulomb  double-counting  correction9 
is  a  similar  contribution  of  —7.24  eV.  The  Ewald 
and  exchange-correlation  correction  energy  differ¬ 
ences  are  13.10  and  0.72  eV,  giving  the  final  energy 
change  (gain)  from  relaxation  of  —1.12  eV.  Thus 
the  eigenvalue  sum  contribution  has  the  same  sign 
as  the  total  energy  difference  in  this  case,  but  other 
contributions  compete  in  importance  with  the  eigen¬ 
value  sum. 

C.  Charge  density 

The  overall  appearance  of  the  charge  density  of 
the  relaxed  HSV  in  the  ( 1 10)  plane  is  similar  to  that 
of  the  ideal  HSV  in  I  and  is  not  shown.  The  differ¬ 
ence  plot  in  Fig.  5  of  p(configuration  3)  —p( ideal) 
indicates  that  the  charge  distortion  accompanying 
atomic  relaxation  is  dominated  by  a  removal  of 
charge  from  between  the  tetrahedron  of  H  atoms 
and  the  addition  of  charge  to  Si  — H  bonding  region. 
This  rearrangement  of  charge  can  be  accounted  for 
qualitatively  by  the  rigid  motion  of  atomiclike  H 
and  Si  charge  densities.  In  addition,  there  is  an  in¬ 
crease  in  charge  in  the  H— Si  backbonding  position 
(more  than  l  a.u.  from  the  Si  nucleus)  which  is  con¬ 
sistent  with  this  ngid-atom  picture.  The  increase  in 
charge  in  the  Si  — Si  bond  likewise  can  be  ascribed 
to  the  2.5%  compression  of  this  bond;  its  asym¬ 
metric  -onn,  however,  suggests  a  distortion  of  the 
bond  charge  due  to  bond  rotation.  The  atomic  re¬ 
laxation  also  results  in  a  decrease  of  the  already 
small  charge  density  in  the  interstitial  region. 


CHARGE  DENSITY  DIFFERENCE 


(HD)  PL-'.NE  ♦  S, 

•  H 


FIG.  5.  Plot  in  the  (110)  plane  of  the  difference  in 
charge  density  between  configuration  3  and  configura¬ 
tion  0;  i.e.,  "relaxed"  minus  "ideal."  Contours  are 
drawn  at  0  iheavy  contour),  +5,  ±  10,  ^20,  and  -50 
electrons/cell. 
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IV.  DISCUSSION 

The  preceding  section  has  described  the  sym¬ 
metric,  constant  volume  relaxation  of  H  and  neigh¬ 
boring  Si  atoms  in  the  HSV.  The  restriction  to  con¬ 
stant  volume  has  compressed  the  H  — Si  and  Si— Si 
bonds  by  7%  and  2.5%  compared  to  their  “equili¬ 
brium”  values,  indicating  the  “relaxed”  HSV  as 
described  here  is  under  pressure.  The  compressibili¬ 
ty  of  crystalline  silicon  and  the  Si  — Si  bond 
compression  (corresponding  to  AK/F=7.5%)  im¬ 
plies  a  local  pressure  of  the  order  of  75  kbar.  A  full 
relaxation  of  the  HSV  (within  the  breathing-mode 
constraint  assumed  here)  would  require  minimizing 
the  energy  with  respect  to  lattice  constant,  then 
again  allowing  internal  atomic  relaxation,  and 
iterating  this  procedure  to  convergence.  It  is  not 
clear  that  the  result  of  such  a  calculation  would  jus¬ 
tify  the  expense,  so  “full  relaxation”  in  this  sense 
has  not  been  attempted. 

As  a  first  approximation  the  “fully  relaxed”  HSV 
(in  the  sense  described  above)  can  be  approximated 
by  the  present  relaxed  atomic  configuration,  but 
with  lattice  constant  expanded  by  2.5%  to  restore 
the  Si— Si  bond  to  its  length  in  crystalline  silicon. 
This  leaves  the  H— Si  bond  compressed  by  4.5%  by 
the  H-H  interactions. 

One  important  property  which  can  *v  sensitive  to 
the  variation  of  lattice  constant  and  to  internal 
stresses  is  the  vibrational  frequency  of  the  system. 
From  the  potential-energy  expression  (4)  and  atomic 
masses  the  breathing-mode  frequencies  are  found  to 
be  3300  and  425  cm-1  for  “optic”  and  "acoustic” 
modes,  respectively.  [These  modes  however  are  not 
normal  modes  of  the  present  supercell  lattice  since 
the  Si(2)  atoms  are  fixed.]  The  high-frequency  H-Si 
stretching  mode  is  50%  above  the  2100-cm-1  band 
assigned  to  H-Si  stretching  modes14-16  in  a-SiH,. 
Volume  relaxation  would  lower  the  calculated  fre¬ 
quency,  perhaps  drastically,  considering  the  large 
but  canceling  electronic  and  ionic  forces  on  H 
(Table  II)  at  equilibrium. 

Experience  in  electronic  structure  calculations  in¬ 
dicates,  however,  that  the  qualitative  features  of  the 
electronic  structure  of  the  HSV  should  not  vary  sig¬ 
nificantly  with  a  2.5%  increase  in  lattice  constant. 
Thus  the  changes  in  the  LDOS  near  the  gap  should 
be  taken  seriously  and  their  implications  given  con¬ 
sideration. 

Brodsky4  has  suggested  a  “quantum-well  model" 
for  a-SiHj,  in  which  it  is  assumed  that  H  — Si  bonds 
deplete  the  valence-band  DOS's  near  the  gap  and  re¬ 
sult  in  localization  of  the  wave  functions  within  0.6 


eV  of  the  valence-band  edge.  Support  for  this 
model  arose  from  the  isolated  HSV  calculations  of 
Di Vincenzo  et  al?  (and  also  from  I),  and  the  model 
could  explain  qualitatively  certain  optical  and  trans¬ 
port  processes.  On  the  contrary,  the  present  study 
shows  that  most  of  the  H  — Si  bonding  states  are  re¬ 
stored  to  the  region  —0.6  <E  <0  eV  by  symmetric 
relaxation. 

The  CPA  HSV  alloy  model  of  Papaconstanto- 
poulos  and  Economou1  (PE)  also  indicates  a  strong 
depletion  of  valence-band  -states  just  below  the  gap. 
Moreover,  this  model  explains  in  a  quantitative 
manner  sl  number  of  properties  of  a-SiHx:  (1)  re¬ 
moval  of  states  from,  and  the  widening  of,  the  gap 
upon  hydrogenation,1  (2)  existence  of  Si  — H  anti¬ 
bonding  resonance  states’  5  just  above  the  gap,  (3) 
the  absorption  coefficient6  in  the  range  2  <  hv  <  3 
eV,  and  (4)  the  position  of  H-related  bonding  states 
in  the  valence  band.1  From  the  “H-site”  and 
“H— Si  bond"  LDOS  in  Fig.  4  it  is  clear  that  sym¬ 
metric  breathing  relaxation  of  the  HSV  removes  pre¬ 
cisely  those  features  of  the  H-related  LDOS  in  the 
CPA  model  which  are  responsible  for  giving  an  inter¬ 
pretation  of  properties  (l)-(i)  above.  Property 
(4) — the  position  of  H-related  photoemission 
peaks — is  less  sensitive  to  the  relaxation  calculated 
here,  with  relaxation  causing  peaks  to  shift  away 
somewhat  from  those  measured  by  von  Roedem 
et  alf  but  more  toward  agreement  with  those 
found  by  Smith  and  Strongin.18 

As  a  consequence,  it  appears  that  the  ideal  HSV 
provides  a  good  model  of  a  local  defect  upon  which 
an  alloy  model  of  a-SiH,  can  be  built,  while  our  re¬ 
laxed  HSV  fails.  This  result,  in  itself,  implies  no 
contradiction,  as  PE  have  emphasized  that  it  is  the 
local  H-Si  chemical  environment  rather  than  the 
geometrical  arrangement  of  the  cluster  of  four  H 
atoms  which  is  the  determining  factor  in  their 
model. 

This  result  may  however  be  interpreted  as  sug¬ 
gesting  only  that  the  present  model  of  relaxed  HSV 
does  not  actually  occur,  but  rather  is  dynamically 
unstable  towards  a  lower-symmetry,  lower-energy 
configuration.  The  present  study  guarantees  only 
that  the  relaxed  configuration  is  stable  with  respect 
to  breathing-type  distortions  with  full  tetrahedral 
symmetry,  and  that  it  is  not  Jahn-Teller  unstable  to¬ 
ward  a  lower-symmetrv  configuration.  However, 
the  restricted  minimum-energy  configuration  found 
here  may  well  be  a  saddle  point  in  the  full  configu¬ 
ration  space  of  atomic  distortions. 

In  such  a  case  the  H  (and  Si)  atoms  are  dynami¬ 
cally  unstable  with  respect  to  motion  perpendicular 
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to  the  (1 1 1)  H  — Si  bonding  directions.  The  fully  re¬ 
laxed  configuration  could  be  a  low-symmetry  one 
which  allows  lengthening  of  the  compressed  H— Si 
and  Si— Si  bonds  toward  more  typical  lengths  as 
well  as  the  avoidance  of  strong  H-H  interactions,  at 
the  cost  of  relatively  low-energy  distortions  of 
H— Si— Si  and  Si— Si— Si  bond  angles.  Avoiding 
the  compression  of  the  H— Si  bond  should  leave  the 
H  LDOS  depleted  in  the  upper  valence-band  region, 
similar  to  that  of  the  ideal  HSV  in  I,  and  therefore 
retain  the  single  feature  (the  local  H-Si  chemical  en¬ 
vironment  discussed  by  PE)  most  essential  to 
modeling  a-SiHx. 

Assuming  a  complete  loss  of  symmetry,  the  con¬ 
figuration  space  consists  of  the  coordinates  of  each 
of  the  four  H  and  four  Si  atoms,  a  total  of  24  vari¬ 
ables.  The  present  approach  of  calculating  forces 
rather  than  energies  provides  the  only  reasonable 
approach  to  finding  the  lowest-energy  configuration 
in  such  a  case.  However,  the  loss  of  symmetry 
leads  to  a  large  increase  in  computational  effort 
(from  k-point  sampling),  and  nothing  is  known  at 


present  about  the  number  of  configurational  itera¬ 
tions  to  be  expected  to  reach  equilibrium  in  such  a 
large  parameter  space.  Nevertheless  restricted  re¬ 
laxations  along  these  lines  may  be  studied  in  the 
near  future. 
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CALCULATIONS  OF  TRANSPORT  PROPERTIES  IN  a-SiiH 

W.E.  Pickett,  D.A.  Papaconstantopoulos  and  E.N.  Econonou 
Haval  Hesetu'ch  l.aboracory,  Washington,  1C  2021 S,  1.5. A. 


Abs  t  rac  t .  -  Ue  have  used  the  coherent  potential  approximation 
to  calculate  the  d.c.  conductivity  for  a  mode!  of  hydrogenated 
amorphous  silicon.  The  results  show  that  the  mobility  of 
electrons  is  greater  than  that  of  holes.  This  appears  to 
be  related  to  the  strong  hydrogen  component  of  the  density 
of  states  just  above  the  gap. 

The  incorpo ration  of  hydrogen  into  amorphous  Si  has  been  found 
to  change  its  electronic  structure  and  transport  properties  drastic¬ 
ally.  The  density  of  states  in  the  band  gap  is  reduced  by  several 
orders  of  magnitude  (1),  thereby  allowing  the  "adjustment",  by  doping 
(2)  with  impurities,  of  the  Fermi  level  Ep  which  is  otherwise  pinnea 

near  the  center  of  the  gap  in  a-Si.  Spectral  features  related  to 
hydrogenation,  such  as  the  removal  of  states  (1)  from  the  jap.  the 
widening  of  the  optical  gap  (1,3),  and  the  alteration  of  both  the 
valence  (4)  and  conduction  densities  of  states,  have  been  studied  by 
a  variety  of  theoretical  methods  (5-8). 

The  model  proposed  by  two  of  the  authors  (6,7)  accounts  for  the 
data  mentioned  above  as  well  as  showing  substantial  agreement  (7) 
with  the  optical  absorption  data  of  Cody  el  aj_  (2).  In  this  paner  we 
begin  a  theoretical  study  of  the  transport  properties  of  a  -  5 • : -  by 
calculating  the  d.c.  conductivity  for  this  model  system. 

The  model  consists  of  an  effective  lattice  whose  sites  have  a 
probability  c  of  being  vacant  and  probability  1-c  of  being  occupied 
by  a  Si  atom.  This  Si-vacancy  alloy  provides  a  model  which  exhibits 
many  of  the  properties  (6,7)  of  a  - S 1 .  In  addition,  hydrogen  atoms 
(one,  two,  three,  or  four)  are  incorporated  around  the  vacant  site 
such  as  to  bond  with  neighboring  Si  dangl'ing  bonds,  and  the  resulting 
spectral  properties  reproduce  many  trends  evident  in  the  experimental 
data  on  hydrogenation  of  a-Si.  The  model  is  based  on  a  S 1  a t e r-Ko s t er 
(SK)  Hamiltonian  fit  (3)  to  the  pseudopotential  band  structure  of  Si, 
and  the  Coherent  Potential  Approximation  (CPA)  is  invoked  to  t'-eat 
the  disorder.  The  c a  1 cu 1  a t i o na 1  details  have  been  given  e'sewne^e 
(6,7). 


In  the  present  work  we  treat  only  the  fully  hydrogenated  case  of 
four  hyorogens  per  vacant  Si  site.  The  resulting  system  shows  a  very 
low  density  of  states  in  the  band  gap,  sinilar  to  data  on  a-Si:'-. 
To  evaluate  the  conductivity  o  we  begin  with  the  zero  frequency  kun  o- 
Greenwood  formula 


2  e2fl 

*2nm2 


/dE  (■  tf) 


T  r  <  px  ImG  ( E  +  i  ; 


) p  ImG ( E ♦ i i / > . 


where  <...>  denotes  the  configuration  average  and  f  is  tne  : 
function.  The  evaluation  of  this  expression  in  a  manner  consistent 
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with  the  CPA  has  been  discussed  in  detail  elsewhere  (9).  By  ignoring 
the  vertex  corrections  arising  from  correlations  in  intermediate 
scatterings,  the  Green's  functions  G  can  be  averaged  separately, 
resulting  in  the  effective  medium  function  G.  The  trace  operator 
Tr  in  Eq.  (1)  includes  a  Sum  over  the  Brillouin  zone  as  well  as  a 
matrix  trace  over  the  basis  states,  which  we  have  taken  as  the  Bloch 
sum  of  s-  and  p-orbitals  which  provide  the  representation  for  the  SK 
Hamiltonian  (8)  of  silicon.  In  this  basis  the  momentum  (px)  matrix 
elements  are  given  by 

(  px/m )  i  j  *  'fi'1  [H(k+6kx)-H(k)]u/6kx  (2) 


which  we  have  evaluated  numerically  by  using  a  finite  wavevector 
displacement  ok  .  If  a  transformation  is  made  to  band  states  wh-ich 
diagonalize  H  §nd  give  eigenvalues  Efc,  the  diagonal  elements  of 

Px/m  become  the  band  velocity  d  E  ^  /  d  (If  k  )  .  However,  preliminary 
evaluation  of  c  retaining  only  band 'diagonal  components  of  Px/m  has 

shown  that  the  full  matrix  nature  of  this  operator  must  be  retained. 
As  opposed  to  previous  CPA  evaluations  (10)  of  o  which  relied  on 
single  band  models,  we  emphasize  that  the  present  calculation  retains 
the  full  multiband  character  of  our  model  of  a  - S  i  :  H . 


The  zone  integration  of  our  expression  for  o(E), 


o(E)  = 

t  On l 


i  jinn 


1  i  j  (  IrnGjm)P 


xWIfflGni> 


(3) 


was  done  by  an  analytic  tetrahedron  approach  (11).  The  irreducible 
zone  (1/4 St h )  was  divided  into  4000  equal  volume  tetrahedra  (compris¬ 
ing  916  distinct  vertices),  with  the  integrand  being  assumed  linear 
within  each  tetrahedron.  The  integration  within  each  tetrahedron  can 
be  evaluated  analytical  ly  ,  and  the  result  is  summed  over  all  4000 
tetrahedra.  The  same  mesh  was  used  to  solve  for  the  CP  A  Green's 
function  U  and  the  related  density  of  states  N(E).  At  zero  tempera¬ 
ture  the  thermal  redistribution  (-3f/9E)  reduces  to  a  5-function 
located  at  Ep,  and  Eq.  (1)  reduces  to  Eq.  (3). 

The  results  for  N  ( E )  and  o(E)  near  the  gap  region  are  shown  in 
Fig.  1.  We  find  that  states  very  near  the  top  of  the  valence  band  at 
E  -  -0.5  eV)  are  only  poorly  conducting,  as  the  initial  increase  in 
3(E)  for  negative  energies  occurs  belowthat  of  N(E).  At  the  top  of 
the  gap,  however,  both  o(E)  and  N(E)  show  the  same  "onset"  at  -  0.S5 
eV.  (The  energy  zero  has  been  set  at  the  bottom  of  the  gap  of  the 
reference  crystalline  Si  system  (8).) 


The  comparison  of  the  contribution  of  electron  and  hole  spates 
to  conductivity  can  be  made  if  we  define 


J  (E )  ~  const.  x<px(E)>N^(E)  (4) 

which  is  suggested  by  the  relation  N(E)  ~  I n G ( E  + i :  ) .  The  mean 

2 

square  matrix  element  <PX(E)>,  which  describes  everything  except 

density  of  states  effects,  is  5  to  15  times  larger  for  electron 
states  than  for  hole  states.  This  difference  reflects  the  character 
of  the  states  which  are  involved:  hole  states  derived  from  bonding 

combinations  of  Si-Si  and  Si-H  orbitals,  and  electron  states  derived 
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Energy  (eV) 

Fig.  1  The  zero  temperature  conductivity  a ( E )  {full  line)  and  density 
of  states  11(E)  (dashed  line)  near  the  gap  region  for  the  present 
model  of  fully  hydrogenated  amorphous  silicon. 


from  anti-bonding  combinations  of  these  orbitals.  Both  Si-H  combi¬ 
nations  are  pulled  down  in  energy  relative  to  their  Si-Si  counter¬ 
parts,  and  this  is  reflected  in  the  ratio  NH(E)/Nt(E)  of  H  local 

density  of  states  N„  to  the  total  Nt  being  three  times  larger  for 
electron  states  than  for  hole  states  {-7*  to  2.5%). 


The  finite  temperature  conductivity  (ignoring  thermal  lattice 
disorder)  can  be  calculated  from  the  relation 


o(T) 


3f(E-EF) 

___ 


(5) 


The  results  are  shown  in  rig.  2  for  various  values  of  Ec  within  the 
gap,  c orrespond i ng  physically  to  different  doping  levels.'  The  calcu¬ 
lated  values  of  log  a  vs.  1/T  reproduce  the  approximately  linear 
behavior  found  in  intrinsic  and  doped  samples  of  a-Si:H  (2,12). 


Fig.  2  -Logarithm  of  a  ( ull 
cm)"1  vs  inverse 
temperature  for 
various  placements  of 
Ep,  corresponding  to 
different  doping 
levels.  Solid  lines 
denote  values  of  Ep 

in  the  lower  half  of 
the  gap  (-0.4,  -0.2, 
0.0  eV ) ,  while  d ashed 
lines  are  used  for 
values  of  E-  in  the 
upper  part  of  the  gap 
(0.2,  0.4,  0.6  e  V)  . 


103/T  (K-1) 
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We  present  the  first  application  of  the  c oh e ren t- pot en t i a  1  - 
approximation  In  the  evaluation  of  optical  absorption  a(E)  ot  a- S i H x 
using  a  realistic  multi-band  model.  The  optical  gap  is  larger  than 
the  calculated  density  of  states  gap  and  the  theoretical  a(E)  agrees 
well  with  experimental  data.  These  results  suggest  that  <x(E)  is 
determined  primarily  by  the  local  H-Si  configuration  and  short  range 
order,  but  that  it  is  insensitive  to  the  particular  long  range  order, 
which  is  not  included  in  our  model. 
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In  spite  of  several  recent  developments  in1 'the  theoretical 
understanding  of  the  electronic  spectrum  of  a-SiHx,  a  gap  remains 
in  the  interpretation  of  measured  quantities'”^  such  as  the  optical 
absorption  coefficient  a,  photoconductivity,  photoluminescence,  etc. 
These  properties  depend  not  only  on  the  electronic  density  of  states 
(DOS),  H ( E ) ,  but  also  upon  transition  matrix  elements  involving  the 
eigenfunctions  of  the  system.  Since  real i s t i c  matrix  elements  are 
much  more  difficult  to  obtain  theoretically  than  is  N(E),  the  usual 
approach  is  to  assume  that  the  matrix  elements  are  constant,  i.e. 
energy-independent.  With  this  assumption  a(w)  is  proportional  to 
the  convoluted  DOS  (CDOS) 

ft 

a (u>  )  -ui”1  I  Nv(E)Nc(E+lMdE  (1) 

Ef--Kw 

(v,c  denote  valence,  conduction  band)  and  the  resulting  optical  gap 
E°pt coi nc i des  with  the  DOS  gap  E®®^.  In  this  paper  we  show  that 
both  the  magnitude  of  a  and  the  measured1’^  E°pt  of  fully  hydrogenated 
a-SiHx  are  reproduced  by  our  calculations  if  we  remove  the  assumption 
of  constant  matrix  elements,  and  a  better  u n de rs t a nd i ng  of  the 
optical  process  is  thereby  obtained. 

In  the  present  work  we  calculate  a(w)  within  a  recently  devel- 
4  5 

oped  multiband  model  ’  of  a-SiHx>  These  calculations  are  the  first 
application  of  the  coherent-potent  i  a  1  -a  pprox  imat  i  on  (Cf’A)  in  the 
evaluation  of  optical  properties  of  a  real  material  using  a  realistic 
multiband  system;  previous  calculations®  of  the  optical  conductivity 
have  been  limited  to  simple  single  or  double  band  models.  The 
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significant  aspects  of  our  work  include:  (i)  wftJhout  adjustable 

parameters  we  obtain  good  agreement  with  experiment*  for  a(w),  both 

in  energy  dependence  and  absolute  magnitude,  and  (ii)  we  find  that 
cOpt  v  cDOS 


>  E 


is  due  primarily  to  remanent  conservation  of  the  magni- 


9  9 

tude  of  the  wavevector  resulting  from  short-range  order. 

The  optical  absorption  coefficient  is  given  in  terms  of  the 

optical  conductivity  o(u) )  and  index  of  refraction  n(u>)  by 

a(iu)  =  4uRe  0  (o))/cn(aj).  (2) 

where  c  is  the  speed  of  light.  The  weak  u -dependence  of  n  has  been 
1 


and  is  included  in  the  calculations  below. 
7 


measured 
use  the  general  expression 

E 


For  o  (  u) )  we 


Rea(u)  =  2~l .  5  f  dE  Tr<pImG(k,E+iO)pImG(t ,  E+fiw+ i 0 ) >  , 

ufinrw  J 


E  p-lfw 


(3) 


2e‘ 


2 

irflm  a) 


F 

w  dE  Tr|px<ImG(k,E  +  iO)>px<ImG(k  ,  E^t  iG  )  >|. 


E  p-fiio 


where  Q  is  the  normalization  volume,  px  is  the  x-component  of  the 
momentum  operator  and  G  is  the  Green's  function  of  the  system.  The 
brackets  <....>  indicate  an  ensemble  average  over  the  disorder  and  A 
denotes  that  A  is  an  8x8  matrix.  For  G  we  use  CPA  results  for  a 

4 

highly  successful  model  a-SiHx  system  discussed  elsewhere.  The 
model  consists  of  a  tight-binding  description  of  20  at.X  H  bonded  to 
5t  Si  vacancies  in  random  positions  on  a  Si  lattice.  The  second 


Q 

approximate  relation  (3)  results  from  ignoring  vertex  corrections 
to  the  exact  expression,  and  allows  the  evaluation  of  Reo  directly 
from  the  CPA  Green's  function  <G>.  The  momentum  matrix  is  given  by 


dHQcl 

dk 


(4) 


where  H  is  the  virtual  crystal  tight-binding  Hamiltonian. 

The  results  for  Reo(u-o.Ep)  have  been  presented  previously^ 
as  a  function  of  Ep.  There  it  was  found  that  contributions  from 
off-diagonal  matrix  elements  of  px  and  <G>  were  negligible  (typically 
~2%) ,  consistent  with  an  " 1 ntraband"  picture  of  the  d.c.  conductivity. 
In  crystals  absorption  at  finite  frequency  occurs  due  only  to  inter¬ 
band  transitions,  but  disorder  admits  the  possibility  of  " i ntraband" 
transitions.  Whereas  it  will  be  seen  that  the  k- dependence  of  the 
matrix  elements  is  necessary  to  understand  the  energy  gaps,  the  full 
matrix  summations  in  the  trace  in  Eq.  (3)  are  necessary  to  obtain  the 
correct  magnitude  of  a. 

Our  result  for  a(w),  calculated  from  Eqs.  (2-4)  at  T=0  (no 
thermal  disorder)  with  no  adjustable  parameters  is  compared  in 
Fig.  1  (as  the  dashed  line)  in  the  form  (aE)1^  vs  E  (E=1W)  with  the 
experimental  data  of  Cody  £t  aj,.  *  (symbols)  and  with  the  COOS  expres¬ 
sion  of  Eq.  (1)  (dotted  line)  corresponding  to  energy- 1 ndependent , 
non-momentum  conserving  matrix  elements  of  px.  In  the  range  2.1  ± 


E  <  3.0  eV,  (a£) 


1/2 


is  in  excellent  agreement  with  the  data  apart 


from  an  overall  factor  of  1.5.  In  the  region  1.6  <.  E  <.  2.1  eV  the 
calculation  leads  to  higher  absorption  (relative  to  the  higher  energy 
region)  than  is  found  for  the  sample  shown  in  Fig.  1.  This  could 


be  due  to  the  lack  of  full  atomic  disorder  or  td'the  absence  of 
off-diagonal  disorder,  either  of  which  would  require  a  more  sophis¬ 
ticated  CPA  model. 

Another  possibility,  which  we  now  address,  is  that  the  approxi¬ 
mation  of  Eq.  (3)  of  neglecting  vertex  corrections  is  breaking  down 
in  this  region.  Examination  of  contributions  to  a(w)  for  fiu  <  2eV 
from  various  regions  of  k-space  shows  that  this  absorption  is  due 
entirely  to  transitions  for  which  either  (k,E)  or  (l.E+'fiw)  is  rather 
far  from  a  broadened  "band"  of  a-SiHx.  In  terms  of  the  spectral 
density  S  given  by 

S(k,E)  =  -  w ” 1 T rG(k,E+iO)  ,  (5) 


either  initial  or  final  state  of  the  transition  lies  where  S  is 
small.  The  CPA,  however,  is  designed  primarily  to  provide  average, 
macroscopic  properties  of  the  system,  and  it  may  be  overestimating 
S ( k , E )  in  regions  far  away  from  bands  which,  for  example,  contribute 
little  to  N(E).  In  a  region  where  the  influence  of  a  single  band  m 
with  energy  E-j^m  is  dominant,  S(k,E)  assumes  the  form 


sfk.E) 


I  I Imrm(E)l 


(6) 


since  for  the  tight-binding  CPA  we  find  the  self-energy  E  (k,E)  in 

the  band  representation  is  nearly  k-i ndependent.  Therefore  S(k,E)  is 

lorentzian  in  E^m  at  fixed  E.  Observation  of  first  principles  muffin- 

q 

tin  CPA  spectral  functions  however  suggests  a  much  faster  than 
Lorentzian  fall-off  of  S(k,E);  i.e.  {Imrm(k,E)(  should  dec reas e  as  k 
leaves  the  vicinity  of  band  m. 
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Following  this  line  of  reasoning  we  modified  the  CPA  self¬ 
energy  to  conform  to  this  more  physical  behavior  by  defining  a  k- 

~  10 

dependent  on-site  self-energy  z  by 


Imrs<p(k.E) 


's*p  (E-E£n)2  +  r2 


ImIs,p(E> 


(7) 


r  =  c  |lmrs(E)+  3 Imrp( E ) I  /4,  (8) 

-► 

where  Ekn  is  the  eigenvalue  at  k  nearest  to  E.  The  real  part  of 
I  Is  relatively  less  important  in  the  calculation  and  was  left  un¬ 
changed.  The  constants  c$  and  cp  were  determined  iteratively  by 
the  requirement  that  the  derived  "modified  CPA"  (MCPA)  DOS  N$  ,  Np 

differ  by  less  than  10%  from  the  CPA  values.  The  value  c  =  0.2  was 
used,  but  with  this  self-consistency  constraint  on  N$  and  Np  the 

results  are  weakly  dependent  on  the  choice  of  c. 

The  effect  upon  S(k,E)  of  this  modification  is  shown  in  Fig.  2 

for  k  along  the  a  direction  and  for  two  energies,  -0.8  eV  and  1.1  eV, 

which  contribute  to  absorption  at  'flu)  *  1.9  eV.  S  peaks  near  the 
eigenvalues  of  the  underlying  crystalline  Hamiltonian.  Far  from  a 
peak,  S  in  the  MCPA  is  typically  reduced  by  an  order  of  magnitude, 
which  in  turn  reduces  the  low  frequency  absorption.  By  constructing 
this  more  physical  form  for  <ImG>  by  our  MCPA  procedure,  we  obtain 
corrections  to  a  which  must  be  included  in  the  vertex  corrections 

neglected  in  Eq.  (3).  The  result  is  shown  in  Fig.  1  as  the  full 

1/2 

line.  The  energy  dependence  of  (aE)  calculated  from  the  MCPA  is 
slightly  steeper  at  higher  energy  than  the  experimental  values  of 
Cody  et  aK ,  but  the  behavior  below  2  eV  is  much  improved  over  the 
CPA  results. 
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Comparison  of  the  calculated  value  of  (aEj^^'-'to  that  arising 
from  a  simple  CDOS^  demonstrates  that  E ° ^ ^  >  E^^  (Fig.  1).  Linear 
extrapolation  of  (a£)  '  to  zero  from  the  linear  high  frequency 
region  leads  to  a  value  E°pt  as  1.8  eV  (CPA)  or  2.1  eV  (MCPA),  whereas 
ED0S  =  1.4  eV  from  the  COOS.  Cody  et  al.1  have  shown  that  E0^* 

g  —  —  g 

varies  from  1.5  to  2.0  eV  with  increasing  hydrogen  content.  It  is 
encouraging  that  our  model,  which  contains  20X  H,  gives  an  optical 
gap  of  the  order  of  2  eV,  in  excellent  agreement  with  data1  on 
samples  with  similiar  H  content.  By  making  the  replacement  px  — > 
diagonal  constant  in  Eq.  (3)  we  obtain  a  similiar  frequency  depend¬ 
ence,  indicating  the  widening  of  the  optical  gap  is  due  primarily  to 
► 

k-space  correlations  of  the  electron  and  hole  in  the  excitation 
process  rather  than  an  explicit  energy  dependence  of  px-  Physically 
this  remanent  k  conservation  is  due  to  the  short-range  order  in 
a-SIH  .  Note  that  our  value  E°^*  s#  2  eV  is  much  smaller  than  that 

*  9 

4 

for  crystalline  Si  ('3  eV)  in  spite  of  the  widening  of  the  gap  due 
to  lowering  of  the  valence  bands  upon  incorporation  of  H.  All  of  the 
absorption  below  3  eV  arises  from  the  disorder  in  the  model,  and  the 
close  agreement  with  experiment  indicates  both  that  the  underlying 
"crystallinity"  of  our  model  is  not  a  serious  defect  and  that  knowl¬ 
edge  of  the  long  range  disorder  In  a-S1Hx  is  not  essential  for 
describing  a. 

We  performed  two  additional  calculations  which  illustrated  other 
features  of  our  model.  First,  we  reduced  the  density  of  fully  hydro¬ 
genated  vacancies  from  5%  to  2.5%.  This  Introduced  two  competing 
effects;  the  decrease  of  the  H/S1  ratio  reduced  Eg^$,  but  since  the 


disorder  was  also  severely  reduced,  E°pt  actual ly • wi dened  somewhat, 
demonstrating  the  approach  to  the  crystalline  Si  limit  E°pt  a-  3  eV. 
The  net  effect  was  to  reduce  (aE)*^  by  a  factor  of  two  below  E  * 
2.6  eV.  Secondly,  we  replaced  1/4  of  the  H  atoms  in  the  original 
calculation  by  dangling  bonds.  This  caused  a  dramatic  increase  of 
(aE)*^  below  2  eV  due  to  both  the  increased  disorder  and  dangling 
bond  gap  states,  and  resulted  in  a  much  smaller  E°pt.  Although  this 
result  follows  the  trend  in  the  data  when  the  H  content  is  reduced. 
It  should  not  be  expected  to  agree  numerically  with  the  experimental 
E°pt- vers us-H-content  line*  because  the  dangling  bonds  have  not  been 
allowed  to  reconstruct.  "Fully  hydrogenated"  a-SiHx  does  not  recon¬ 
struct  significantly  because  all  potential  dangling  bonds  have  bonded 
with  hydrogen.  Therefore  our  model,  when  fully  hydrogenated,  corre¬ 
sponds  semiquant i tat i vely  to  fully  hydrogenated  samples  (as  long  as 
localized  band  tail  and  gap  states  are  not  crucial),  but  it  does  not 
have  the  flexibility  to  model  situations  with  reduced  H  content  to 
the  same  accuracy. 

To  summarize,  we  have  presented  here  realistic  microscopic 
calculations  of  the  absorption  coefficient  of  a  model  for  a-SiHx. 
The  calculated  absorption  coefficient,  obtained  without  adjustable 
parameters,  is  within  a  factor  of  two  of  the  experimental  value, 
which  is  Itself  sample  dependent.  Considered  alongside  the  other 

4 

successes  of  the  model,  this  agreement  indicates  that  knowledge 
of  the  precise  atomic  positions  are  not  necessary  for  the  understand¬ 
ing  of  many  properties  of  a-SiHx  for  which  band  tail  or  gap  states 
are  not  intimately  involved. 


(eV/W 
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Fig.  1  Comparison  of  calculated  absorption  coefficient  a  with 
measured  value  (data  points)  of  Cody  et^  aj_.  ,  Ref.  1,  plotted 
as  (aE)*^  vs  E.  The  CPA  and  MCPA  results  were  multiplied 
by  1.5  and  2.4,  res pect 1 vely ,  to  normalize  to  the  experimental 
data  at  2.9  eV. 
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Fig.  2  Semilog  plot  of  CPA  and  MCPA  spectral  densities  S(k,E)  at  E  * 
-0.8  eV  and  E  «  1.1  eV  for  IT  ||  (100).  The  MCPA  value  is 
roughly  one  order  of  magnitude  smaller  than  the  CPA  value 
far  from  the  peaks. 
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